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PREFACE 


This  report  contains  an  account  of  the  principal  results  of  an  invest- 
tigation  into  new  stoichiometric  compounds  of  Nd  over  a two-year  period. 

The  overall  objective  of  this  program  was  to  develop  new  high  concentration 
Nd  materials  for  miniature  lasers  emitting  at  1.06u.  At  the  present  time 
there  are  only  a few  such  compounds, ’ Nd  pentaphosphate  being  the  most 
widely  known.  Although  remarkable  in  its  efficiency,  relative  lack  of 
concentration  quenching  and  high  gain,  the  growth  of  pentaphosphate  crystals 
has  proved  to  be  difficult.^  The  crystal  growth  problems  of  the  penta- 
phosphate coupled  with  the  possibilities  for  improved  performance  by  new 
materials  were  motivating  factors  in  our  research. 

The  publication  of  a paper’  on  cubic  (ordered  perovskite)  Cs2NaNdClg 
in  1974  gave  us  the  idea  that  the  remarkable  properties  of  this  material 
could  be  used  as  a starting  point  for  the  development  of  a new  class  of 
stoichiometric  compounds  possibly  superior  to  anything  known  at  the  time. 

The  unusual  property  of  Cs2NaNdClg  is  the  exceedingly  long  lifetime  of 
luminescence  (1.2  ms)  in  the  undiluted  compound.  The  length  of  the  decay 
was  ascribed  to  a large  Nd-Nd  separation,  no  sharing  of  nearest  neighbors 
and  the  presence  of  a center  of  inversion  at  the  Nd  site. 


As  reported  in  Ref.  3 the  new  material  is  of  no  particular  interest 
to  miniature  laser  technology  because  of  a low  emission  cross  section 
(and  therefore  gain)  and  poor  stability  in  air.  However,  the  isolating 
properties  of  the  lattice  permitting  such  a long  lifetime  gave  us  an  idea 
of  how  concentration  quenching  can  be  suppressed.  The  initial  course  of 
our  work  was  therefore  directed  at  maintaining  a similar,  ordered  perovskite 
lattice  but  removing  the  inversion  center.  It  was  thought  that  a compound 
with  these  characteristics  might  result  in  a enhancement  of  the  transition 
probability  while  still  maintaining  the  lack  of  concentration  quenching. 


Our  work  in  the  KF-LiF-NdF^  system  did  not  produce  an  ordered  perov- 
skite structure  with  a distorted  site  symmetry.  However,  a new  and  inter- 
esting phase  was  discovered  which  may  prove  a valuable  addition  to  the  list 
of  isolated  Nd  compounds.  Since  the  completion  of  the  contract  period  the 
new  compound,  KgNdLi2Fj^Q,  has  been  isolated  with  properties  that  challenge 
NdPgOj^^  as  a miniature  laser  material.  Prompted  by  a desire  to  eliminate  " 

active  or  passive  Q-switching  while  still  retaining  fast  and  high  peak  power  Section 
pulse  operation,  we  explored  the  dynamics  of  laser  operation  in  a system  ’*‘*®** 

having  two  centers  of  similar  spectral  characteristics  but  having  different 
lifetimes.  A system  of  that  kind  could  be  pictured  as  containing  a low 
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threshold  subsystem  based  on  the  high  gain  centers  and  an  energy  storage  sub- 
system based  on  the  low  gain  centers,  all  in  one  material.  As  it  turned  out 
the  merits  of  this  scheme  are  very  limited.  However,  this  effort  gave  rise 
to  a new  and  comprehensive  analytical  treatment  of  pulse  operation  in 
optically  pumped  lasers. 

This  report  is  organized  as  a series  of  Sections  each  largely  self- 
contained  and  dealing  with  various  phases  of  the  work.  We  summarize  the 
main  achievements  by  presenting  an  abstract  for  each  Section  in  turn: 


SECTION  I 

NEW  LUMINESCENT  Nd^'*'  COMPOUNDS  FOR  LASER  APPLICATIONS; 

I.  ORDERED  CESIUM  SODIUM  NEODYMIUM  HALIDES 

A number  of  polycrystalline  elpasolite  halides  of  the  general  formula 
A2BB'Xg  have  been  prepared  and  investigated  as  possible  candidate  materials 
for  miniature  Nd  lasers.  The  chief  base  materials  have  been  CS2NaNdClg  and 
CS2NaNdFg.  X-ray  diffraction  and  emission  spectroscopy  have  been  used  to 
follow  the  effects  of  various  attempts  to  distort  the  0^  site  symmetry  in 
the  parent-ordered-perovskite  materials.  The  combination  of  excitation 
line  narrowing  (using  monochromatic  excitation)  and  time  resolution  (using 
a boxcar  integrator  gated  at  selected  delay  times)  have  been  very  effective 
in  the  elucidation  of  the  emission  spectra  of  elpasolite  fluorides. 

All  the  fluorides  investigated  contain  two  Nd^''^  sites.  At  the  highest 

Nd^^  concentration,  one  has  a fluorescence  lifetime  of  about  3 us;  the  other 

has  a much  longer  lifetime  of  2 ms.  The  short  lifetime  species  is  badly 

quenched  and  has  been  ascribed  to  NaNd  Y,_  F.  impurity  in  Cs.NaNd  Y,  F,. 

The  quenching  constant,  k^^,  is  more  than  two  orders  of  magnitude  higher 

than  for  stoichiometric  Nd  materials.  The  long  lifetime  centers  show  a 

much  smaller  deviation  of  k (a  factor  of  3 or  4)  from  other  stoichiometric 

cn 

materials.  The  emission  spectrum  of  the  long  lifetime  site  has  been  dis- 
cussed in  terms  of  nearly  symmetry  but  with  the  presence  of  some  lines 
which  are  likely  of  vibronic  origin. 

The  design  of  a stoichiometric  material  based  on  the  elpasolite 
structure  but  with  increased  radiative  probability  is  evidently  much  more 
difficult  than  originally  anticipated.  The  various  attempts  to  reduce  the 
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site  symmetry  by  cationic  and  anionic  substitution  have  been  largely 
unsuccessful.  A further  difficulty  has  been  the  sensitivity  of  cesium 
elpasolites  (both  chlorides  and  fluorides)  to  moisture.  This  work  does 
indicate  that  a better  understanding  of  the  quenching  process  would  be 
helpful  in  future  materials  design.  Finally,  the  usefulness  of  the  methods 
of  investigation  are  shown  both  in  this  section  and  very  successfully  in 
Section  2,  which  covers  the  investigation  of  the  KF-LiF-NdF^  system. 


SECTION  II 

NEW  Nd^'*’  COMPOUNDS  FOR  LASER  APPLICATIONS; 

II.  POTASSIUM  LITHIUM  NEODYMIUM  FLUORIDES 

The  KF-LiF-NdFj  system  was  investigated  with  the  objective  of  deter- 
mining whether  a new  stoichiometric  Nd  material  might  be  found.  With  the 
aid  of  selective  excitation  with  a dye  laser,  the  fluorescence  spectrum  of 
a new  interesting  phase  was  obtained.  The  fluorescence  lifetime  for  Nd  ion 
in  this  phase  was  found  to  be  280  ps,  which  is  more  than  twice  that  of  the 
leading  candidate  material  for  miniature  lasers,  neodymium  pentaphosphate 
(NdP^Oj^^) . A comparison  of  the  emission  spectra  of  single  crystalline 
NdP^Oj^^  and  the  polycrystalline  potassium  neodymium  lithium  fluoride  show 
the  two  to  be  very  similar  with  the  NPP  emission  at  l.OSy  a bit  more  in- 
tense, but  the  KNLF  phase  exhibits  a much  better  branching  ratio.  Prelim- 
inary dilution  studies  with  cerium  ion  indicates  that  concentration  quenching 
may  be  much  less  in  KNLF  than  in  NPP. 

Since  the  end  of  the  contract  period  some  additional  information  has 
been  obtained  in  cooperation  with  S.R.  Chinn  and  H.  Y-P.  Hong  of  the  MIT 
Lincoln  Laboratory  with  the  use  of  single  crystals."*  The  composition  of 
KNLF  has  been  found  as  a consequence  of  an  x-ray  structural  determination 
to  be  K^NdLi2Fj^Q.  The  unit  cell  is  orthorhombic  with  the  space  group  Pnma. 
The  lattice  parameters  are:  a = 2.065  nm,  b = 0.778  nm,  and  c = 0.690  nm. 
Nd^^  is  coordinated  to  eight  fluoride  ions  with  the  point  symmetry  m.  The 
Nd^^  concentration  is  3.61  x 10^^  cm”^  and  the  closest  Nd-Nd  distance  is 
0.673  nm.  Laser  emission  has  been  obtained  with  a platelet  pumped  by  a 
dye  laser.  This  is  the  first  high  concentration  fluoride  material  in  which 
lasing  has  been  observed. 
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SECTION  III 


CONCENTRATION  QUENCHING  IN  Nd^’*'  STOICHIOMETRIC  MATERIALS 

The  basic  experimental  facts  concerning  the  luminescence  of  stoi- 
chiometric compounds  have  been  reviewed  and  further  confirmed  in  the  case 
of  Nd^Laj^_^PgOj^^ . In  particular  we  draw  attention  to  a quenching  rate 
which  is  linear  in  x and  the  exponential  nature  of  the  decay  for  all  x. 

We  point  out  that  the  conventional  models  of  concentration  quenching  based 
on  energy  migration  governed  by  a diffusion  equation  can  only  account  for 
observations  if  it  is  assumed  that  the  stoichiometric  materials  (specifi- 
cally the  pentaphosphate)  contain  a fixed  impurity  concentration  ("vlO^^  cm 
independent  of  x. 

This  model  is  considered  artificial  and  another  one  is  proposed  based 
on  quenching  by  nearby  Nd-Nd  pairs.  This  too  leads  to  difficulties.  A 
conclusion  is  reached  that  there  is  at  present  no  clear  theoretical  picture 
capable  of  explaining  the  results.  Experiments  are  suggested  for  the 
measurement  of  energy  migration  rates  as  a function  of  concentration. 

These  could  give  rise  to  a major  revision  of  conventional  concepts  of 
energy  transfer  and  concentration  quenching. 

SECTION  IV 

APPROXIMATE  ANALYTIC  APPROACH  TO  LASER  PULSE  PHENOMENA 

We  have  suceeded  in  obtaining  approximate  analytic  expressions  that 
describe  a wide  variety  of  laser  pulse  phenomena  that  can  be  handled  using 
the  laser  rate  equations.  Using  these  expressions  we  can  examine  these 
phenomena  without  recourse  to  the  catalogs  of  numerical  computer  data  that 
have  been  the  traditional  approach  to  laser  pulse  phenomena. 

For  laser  materials  with  long  spontaneous  lifetimes  the  Q-switch  mode 
is  to  be  preferred  since  it  gives  the  highest  peak  intensity.  The  enhance- 
ment of  the  peak  intensity  over  cw  is  approximately  equal  to  the  ratio  of 
the  spontaneous  lifetime  to  the  cavity  lifetime. 

As  the  spontaneous  lifetime  approaches  the  cavity  lifetime  the 
preferred  mode  is  cavity  dumping,  if  the  residual  cavity  losses  are  small 
enough.  In  this  case  the  enhancement  over  cw  is  proportional  to  the  ratio 
of  the  cavity  lifetime  with  zero  output  coupling  to  the  cavity  round  trip 
time . 
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Transient  spike  operation  is  useful  in  roughly  the  same  range  as 
Q-switching.  However,  the  enhancement  factor  over  cw  operation  is  low 
relative  to  that  achievable  with  Q-switching.  Transient  spike  operation 
is  useful  if  one  is  interested  in  moderate  pulse  intensities  and  in 
dispensing  with  extra  intracavity  elements  or  if  one  is  principly  interested 
in  obtaining  narrow  pulses.  If  there  are  lasing  centers  of  different  life- 
times present,  transient  spike  operation  is  dominated  by  the  fast  centers 
and  no  enhancement  of  peak  intensity  results  from  adding  long  lifetimes 
centers  for  energy  storage. 

Where  one  is  interested  principally  in  obtaining  laser  pulses  short 
relative  to  those  of  a pump  laser  transient,  pulse  operation  can  be  used 
to  advantage  even  where  the  spontaneous  lifetime  approaches  the  cavity 
lifetime.  With  this  approach  subnanosecond  pulses  can  be  obtained  in  dye 
lasers . 


SECTION  V 

SPECTROSCOPY  OF  A HIGH  Nd  CONCENTRATION  Al-PHOSPHATE  GLASS 

As  an  extension  of  the  work  on  stoichiometric  materials,  we  have  in- 

21  3 

vestigated  a Nd  glass  containing  up  to  2.7  x 10  Nd  ions  per  cm  , or 

2 3 wt  % of  Nd,  and  thus  approaching  the  concentration  in  Nd  pentaphosphate 
21  3 

(which  is  3.9  x 10  Nd/cm  ).  This  work  was  not  part  of  the  contract 
since  it  was  carried  .out  on  internally  provided  funds.  We  include  it  in 
this  report  in  order  to  draw  attention  to  high  concentration  glasses  as 
possible  candidates  for  minilasers.  Contrary  to  widely-held  beliefs, 
glasses  containing  such  high  concentrations  of  Nd  can  still  show  substan- 
tial efficiency.  The  luminescence  decay  times  were  of  the  order  of  40  ps, 
thus  only  a factor  of  three  lower  than  for  the  pentaphosphate. 

Laser  experiments  performed  subsequently  to  the  writing  of  this  report 
showed  that  thresholds  are  between  one  and  two  orders  of  magnitude  higher 
than  for  the  pentaphosphate,  when  both  are  pumped  with  a tunable  dye  laser. ^ 
From  this  it  is  estimated  that  thresholds  for  f lashlamp-pumped  rods  (of  the 
order  of  1-cm  length)  should  be  below  lOJ. 
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1 . INTRODUCTION 


The  purpose  of  the  present  study  is  to  investigate  new  kinds  of 
materials  with  stoichiometric  or  near  stoichiometric  concentrations  of  Nd 
ions.  These  materials  are  excellent  candidates  for  the  development  of 
miniaturized  lasers,  emitting  around  1.06  pm. 

Recently  excellent  review  articles  describing  the  work  done  to  date 
in  several  laboratories  have  been  published  by  Chinn,  Hong  and  Pierce'  and 
by  Danielmeyer . ^ From  these  the  following  criteria  for  the  development 
of  new  stoichiometric  materials  appear  to  emerge:  (a)  the  distance  between 
Nd-Nd  pairs  should  be  large,  (b)  the  Nd  ions  should  not  share  nearest 
neighbors  and  (c)  the  energy  level  scheme  should  he  such  as  to  prevent 
efficient  pair  relaxation.  While  these  criteria  may  not  be  totally  suf- 
ficient and  (c)  may  not  be  easily  controllable,  (a)  and  (b)  are  eminently 
reasonable  in  terms  of  what  is  known  about  concentration  quenching. 

Table  1,  which  summarizes  the  main  properties  of  known  stoichiometric 
materials  is  reproduced  here  from  Ref.  1.  We  wish  to  draw  particular  at ton 
tion  to  the  column  labeled  "Fluorescence  Lifetime.”  It  will  be  noticed 
that  for  all  materials  with  the  exception  of  the  aluminum  borate,  the  life- 
times at  1%  dilution  (x  = 0.01)  are  of  the  order  of  200  ps  to  400  ys,  a 
value  which  appears  to  be  characteristic  for  electric  dipole  transitions  i.n 
Nd^  . The  aluminum  borate  has  a lower  coordination  and  a stronger  de- 
parture from  inversion  symmetry  than  the  rest.  ' 

The  last  entry  in  the  table  represents  a special  case.  The  unusually 
long  lifetime  of  Cs2NaNdClg  is  due  to  the  presence  of  an  inversion  center 
at  the  Nd  site  site  symmetry)  which  makes  electric  dipole  transitions 

symmetry  forbidden.^  We  have  chosen  the  material,  or  rather  materials 
structurally  similar  to  the  hexachloride  as  the  main  subject  of  this  in- 
vestigation. These  materials  (known  as  elpasolites)  have  the  general 
formula  A2BB'Xg.^  The  cation  A is  larger  than  cations  B and  B'  which  are 
nearly  the  same  size  and  occupy  identical  octahedral  sites  in  an  ordered 
manner.  This  serves  to  double  the  cell  parameter  of  the  simple  cubic  ABX^ 

perovskite  structure  and,  importantly,  to  greatly  increase  the  separation 

3+ 

of  Nd  sites. 

Again  referring  to  Table  1,  we  see  that  as  the  concentration  of  Nd 
increases  toward  a purely  stoichiometric  composition,  the  lifetime  shortens 
This  is  a well-known  manifestation  of  concentration  quenching  which  is 
expressible  in  terms  of  the  lifetime  ratio.  In  an  ideal  material,  this 
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ratio  would  be  equal  to  unity,  i.e.,  there  would  be  no  concentration 
quenching.  The  table  indicates  that  by  this  criterion,  there  should  be 
room  for  about  a threefold  improvement.  The  lifetime  ratio  of  the  hexa- 
chloride  does  not  appear  to  be  particularly  low  and,  in  fact,  is  the  second 
highest  of  the  materials  listed.  There  is,  however,  another  consideration 
which  makes  this  structure  of  particular  interest  and  provides  the  main 
motivation  for  the  work  reported  here. 

The  observed  decay  time  of  fluorescence  is  given  by  the  expression: 


where  and  are  the  radiative  and  nonradiative  lifetimes,  respectively. 
Assuming  that  at  1%  dilution  the  observed  lifetime  is  essentially  due  to 
radiative  processes  (i.e.,  = 4100  ys) , we  obtain  from  Ref.  1 for  the 

stoichiometric  (x  = 1)  material  x^  = 1757  ys) . This  is  an  unusually  long 
time  indicating  a very  low  probability  of  nonradiative  deactivation.  It 
may  represent  a consequence  of  the  large  distance  between  Nd  ions  in  this 
structure  (0.77  nm)  and  the  lack  of  any  sharing  of  nearest  neighbors 
between  Nd  pairs. ^ The  structure  of  the  hexachloride  must  thus  have 
excellent  isolating  characteristics  which  prevent  deactivation. 

Suppose  now  that  this  property  of  the  lattice  can  be  maintained  by 
keeping  interionic  distances  approximately  fixed  while  at  the  same  time 
removing  the  center  of  inversion  from  the  Nd  site.  This  should  be  achiev- 
able, for  instance,  by  substitution  of  a different  halogen  (say  a fluorine 
for  a chlorine  in  the  nearest  neighbor  shell) . Such  a substitution  should 
result  mainly  in  a change  of  x^,  leaving  x^^  unchanged.  Suppose  that  x^  can 
be  reduced  tenfold  to  a value  similar  to  that  encountered  in  noncentro- 
symmetric  environments.  Using  Eq.  (1)  , this  should  result  in  x^  ''^300  ys, 
and  thus  a considerable  improvement  in  the  ratio  t^/x^.  The  idea  of 
manipulating  the  immediate  environment  of  Nd  to  enhance  radiative  (electric 
dipole)  transitions  while  retaining  the  highly  isolating  characteristics 
of  the  expanded  perovskite  lattice  is  thus  an  intriguing  one. 

The  limited  success  which  we  report  here  is  due  to  rather  formidable 
chemical  problems  encountered  in  the  preparation  of  mixed  materials.  The 
results  provide,  however,  useful  insights  into  the  mechanisms  of  quenching 
and  for  further  extension  of  the  same  ideas  to  other  systems. 
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2.  EXPERIMENTAL 


2.1  SPECTROSCOPY 

The  experimental  setup  for  spectroscopy  shown  on  Figure  1 enables  us 
to  perform  a variety  of  measurements.  The  sample  is  located  in  a gas  flow 
dewar  whose  temperature  can  be  regulated  anywhere  between  room  and  liquid 
nitrogen  temperature.  For  ordinary  emission  spectra  a 150W  Xe  arc  lamp 
with  a BG18  filter  is  used  for  front  excitation.  Spectra  are  analyzed 
with  a l/2m  Jarrell  Ash  spectrometer,  and  detected  with  a cooled  photo- 
multiplier. In  the  course  of  this  work,  two  different  (PM2)  tubes  of 
vastly  different  responses  were  used.  Earlier  spectra  were  taken  with  an 
ITT,  FW  118  photomultiplier;  later  ones  with  a Varian.  Since  the  Varian 
response  is  very  nearly  flat  in  the  range  800  nm  to  1080  nm,  the  spectra 
obtained  with  it  show  very  nearly  the  correct  intensity  distribution  between 
the  resonance  region  (850  nm  to  950  nm)  and  the  laser  transition  region 
(1020  nm  to  1080  nm) . The  S-1  response  of  the  ITT  tube  tends  to  emphasize 
the  resonance  region  by  a factor  of  about  2.5.  Since  in  this  paper  we  are 
not  concerned  with  branching  ratios,  no  need  for  spectral  correction  of  the 
traces  was  necessary.  In  all  traces  the  number  given  on  the  extreme  left 
identifies  the  original  chart  recording  and  the  conditions  of  experiment. 

All  charts  starting  with  #130  (for  Nd  samples)  and  #E26  (for  Eu  Scimples) 
were  taken  with  the  Varian  tube.  All  lower  numbers  were  taken  with  the 
ITT  tube. 

An  alternative  excitation  source  is  provided  by  a tunable  dye  laser 
of  conventional  design  pumped  by  a nitrogen  laser.  Depending  upon  the 
desired  excitation  wavelength,  rhodamine  6G,  coumarin  120,  or  coumarin  102 
were  used  as  dyes.  The  line  width  of  the  dye  laser  was  typically  0.1  nm 
to  0.2  nm,  the  repetition  rate  60  pps  and  the  tuning  range  typical  of 
rhodamine  6G.  With  monochromatic  excitation  we  can  obtain  "line-narrowed" 
emission  spectra  and  thus  detect  the  presence  of  different  luminescent 
sites.'*  Often,  however,  line  narrowing  is  not  sufficient  to  separate 
spectra  due  to  different  ionic  environments  because  the  absorption  bands 
of  the  different  species  may  still  overlap.  If  the  different  sites  result 
in  luminescence  characterized  by  sufficiently  different  lifetimes,  one  can 
combine  excitation  line  narrowing  with  time  resolution  of  the  detected 
spectra.  This  is  provided  by  a boxcar  integrator  (PAR  162)  whose  gate  can 
be  set  up  for  different  delay  times.  This  technique  has  been  used  exten- 
sively to  separate  the  spectra  of  sites  giving  rise  to  short  (microsecond) 
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and  long  (millisecond)  lifetimes  in  a series  of  materials  described  in 
Section  2.2.  In  general,  therefore,  three  types  of  emission  spectra  were 
taken;  broadband  (Xenon  arc)  excited  to  take  advantage  of  all  pump  bands; 
monochromatically  excited  (dye  laser)  but  time  averaged  by  the  detection 
system  and  finally  monochromatically  excited  and  time  resolved  (boxcar) . 

The  decay  times  can  be  measured  in  two  different  ways.  Fluorescence 
emanating  from  the  back  side  of  the  sample  is  focused  into  a photomultiplier 
through  an  interference  filter  isolating  the  1.06u  group  of  lines.  The 
output  of  the  PM  can  be  displayed  directly  on  an  oscilloscope  or  sampled 
with  the  boxcar.  In  the  case  of  Nd  , this  is  adequate  since  the  lifetime 
of  the  two  main  groups  of  lines  terminating  on  ^19/2  ^^11/2  states  both 

originate  from  the  same  upper  level  ^3/2  which  in  cubic  symmetry  is  not 
split.  In  a distorted  symmetry,  the  splitting  is  usually  quite  small  and 
the  populations  in  the  two  possible  components  are  thermalized.  Alternately 
it  is  possible  to  measure  the  decay  of  a particular  line  by  using  the  spec- 
trometer and  the  PM2  tube.  The  condition  of  excitation  and  the  measured 
lifetime  (where  appropriate)  are  given  on  the  right  side  of  each  emission 
trace. 

2.2  MATERIAL  PREPARATION 

The  goal  of  the  materials  part  of  this  project  has  been  the  preparation 
of  an  ordered  perovskite  neodymium  compound  in  which  the  Nd"^  site  symmetry 
is  reduced  so  as  to  enhance  radiative  transitions.  The  basic  idea  is  to 

distort  the  Nd^^  environment  of  something  like  Cs_NaNdCl,  in  which  Na^  and 

3+  ^ ^ 

Nd  occupy  octahedral  sites  in  an  ordered  manner  as  shown  in  Figure  2.^ 

One  approach  would  be  to  partially  substitute  the  halide  sphere  such  that 

the  anion  neighbors  are  made  up  of  two  different  elements,  such  as  Br  and 

Cl  . Another  approach  would  be  substitution  of  the  cations  in  the  A or 

B sites  in  the  general  A2BB'Xg  formula.  Finally  a solution  to  the  problem 

was  sought  in  widening  the  range  of  elements  such  that  certain  component 

size  restrictions  for  the  perovskite  structure  are  approached.  For  example, 

the  following  tolerance  factor  has  been  used: 


V2- (R3  t Rj^, 

* . 

where  R^,  R^,  and  R^^  refer  to  the  ionic  radii  in  the  perovskite  formula 
ABX^.  For  cubic  oxide  perovskites  0.8  <t  <0.9,  but  a somewhat  wider  range 
has  been  reported  for  the  noncubic  structures  found  mainly  among  oxides.® 
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Figure  2.  Crystal  Structure  of  Cs2NaNdClg 


Our  approach  has  been  to  first  investigate  samples  which  were  expected 

3+ 

to  be  cubic  and  then  attempt  distortion  of  the  Nd  site.  Of  course,  a 
simple  distortion  of  the  unit  cell  as  indicated  by  x-ray  powder  diffraction 
data  does  not  necessarily  imply  a distortion  of  the  Nd^^  site.  In  order  to 
follow  local  changes  in  the  Nd^^  environment,  emission  spectra  of  samples 
were  routinely  taken  as  reported  in  the  Spectroscopy  and  Spectroscopy  Results 
sections.  The  details  of  the  preparations  of  the  various  perovskite  samples 
are  provided  on  the  following  pages. 

CS2NaNdClg 

This  material  is  the  prototype  of  the  ordered  perovskites  containing 
Nd^^  which  are  of  interest  (see  Figure  2) . The  preparation  of  this  material 
and  the  growth  of  single  crystals  of  ordered  Cs2NaMClg  compounds  (where  M 
represents  a number  of  trivalent  metal  cations)  were  first  reported  by 
Morss,  et  al.^  Samples  were  reportedly  prepared  by  evaporation  of  an  HCl 
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solution  of  the  cations  to  dryness  or  by  combination  of  the  anhydrous  halides 
under  vacuum  or  inert  atmosphere  conditions.  These  preparative  techniques 
were  later  used  by  Tofield  and  Weber  to  obtain  polycrystalline  and  single 
crystalline  Cs2NaNdClg  samples  for  their  fluorescence  studies.® 

We  found  the  preparation  of  Cs2NaNdClg  to  be  rather  difficult,  not  so 
much  from  the  air  sensitivity  of  the  product,  but  because  moisture  is  in- 
troduced readily  by  handling  the  reactants  in  air  (especially  NdCl^) . The 
attempted  formation  of  pure  Cs2NaNdClg  by  the  reaction  of  CsCl,  NaCl,  and 
Nd202  in  a stream  of  HCl  has  demonstrated  the  difficulty  of  the  complete 
conversion  of  rare-earth  oxides  to  chlorides  by  this  method. 

The  alternative  method  of  evaporating  an  appropriate  HCl  solution  to 
dryness  has  been  found  to  result  in  impure  samples,  but  the  method  might  be 
acceptable  if  the  evaporation  is  done  slowly,  or  better,  under  reduced 
pressure.  An  improved  procedure,  greatly  reducing  the  water  content  of 
the  final  precipitate  is  achieved  in  the  following  way.®  The  volume  of  an 
HCl  solution  of  component  chlorides  is  reduced  to  the  point  where  Cs2NaNdClg 
just  begins  to  precipitate.  The  sample  is  redissolved  by  the  addition  of 
methanol,  the  volume  is  again  reduced  to  the  point  of  precipitate  formation, 
and  the  process  is  repeated  several  times  with  HCl  being  bubbled  through 
the  solution.  After  several  iterations,  a final  precipitation  is  made  by 
reducing  the  solubility  of  the  product  by  the  addition  of  isopropanol.  This 
is  dried  under  an  HCl  flow  by  gentle  heating  followed  by  somewhat  higher 
temperatures  to  remove  the  last  traces  of  solvent.  Even  this  did  not  prove 
entirely  satisfactory  because  impurities  were  always  shown  to  be  present 
by  x-ray  powder  diffraction  data. 

The  direct,  and  probably  best,  preparation  of  Cs2NaNdClg  involves  the 
reaction  of  the  component  chlorides  in  a stream  of  HCl,  or  in  a sealed 
ampule.  In  either  case  some  precautions  should  be  taken  to  avoid  sample 
contamination.  An  argon  filled  glove  box  and  vacuum  line  techniques  were 
used  during  and  after  sample  preparations  in  order  to  minimize  atmospheric 
contaminants . 

Cesium  chloride  is  somewhat  volatile  at  preparation  temperatures  of 
800°  to  850°C;  thus,  precise  control  of  stoichiometry  is  a bit  difficult  in 
an  HCl  flow  because  of  loss  of  CsCl.  Therefore,  both  polycrystalline  and 
samples  with  single  crystalline  sections  were  prepared  by  heating  dried 
anhydrous  chlorides  in  sealed,  evacuated  quartz  ampules.  Best  results  were 
obtained,  especially  in  crystal  growth  experiments,  after  the  melting  point 
of  Cs2NaNdClg  was  found  to  be  about  715°C  (the  sample  was  observed  to  be 


10 


liquid  at  720°)  rather  than  800°  to  850°C  as  implied  elsewhere.®  Knowing 
this,  samples  were  heated  above  800°C  for  only  a few  hours  (to  allow  complete 
reaction  of  the  components)  and  then  brought  to  about  730°  to  750°  for  the 
crystal-growth  process.  By  this  procedure,  attaclc  of  the  quartz  ampules 
was  visibly  reduced. 

The  crystal-growth  procedure  employed  the  gradient  freeze  method  in  a 
zoned  vertical  furnace  with  gradients  of  about  5°  to  40°C  per  inch.  Samples 
were  equilibrated  25°  to  50°  above  the  melting  point  for  six  to  ten  hours 
prior  to  controlled  cooling  at  2°  to  5°/hour.  The  quartz  sample  containers 
were  6 mm  to  10  mm  I.D.  with  a conical  tip  drawn  at  about  30° to  45°  angles. 
Best  results  were  obtained  when  a constriction  was  provided  about  one-half 
inch  above  the  sample  tip  in  order  to  reduce  crystal  defects  and  twinning. 
The  samples  were  susceptible  to  attack  by  moist  air,  were  rather  soft,  and 
proved  difficult  to  polish.  Chloride  samples  grown  by  this  method  included 
Cs2NaNdClg,  Cs2LiNdClg,  and  Cs2Nag  g^^9Q  2^^6'  showed  clean  x-ray  patterns 

for  the  ordered  perovskite  phase.  Emission  spectra  showed  evidence  that 
Nd^^  is  located  only  in  an  Oh  symmetry  octahedral  environment. 

Cs-NaNd  Y, 

2 X 1-x  6 

Initially,  preparation  of  Cs2NaNdFg  was  attempted  by  substitution  of 
Cs2NaNdClg  by  various  fluorinating  agents  but  without  notable  success.  It 
was  apparent  that  a preparation  from  the  component  fluorides  or  from  ap- 
propriate carbonates  or  oxides  by  the  action  of  HF  in  a Monel  apparatus 
would  be  best.  To  this  end,  the  apparatus  diagrammed  in  Figure  3 was  con- 
structed. The  gas  line  inputs  were  Monel  and  allowed  the  use  of  argon  as 
a flushing  gas  prior  to  heating  the  sample.  A screw-on  Monel  cap  fitted 
with  a Teflon  washer  on  the  upstream  end  of  the  reaction  tube  provided 
access  for  sample  loading  and  withdrawal  as  well  as  cleaning  prior  to  use. 

A nitrogen  flow  was  maintained  over  the  section  of  the  outer  wall  of  the 
Monel  tube  inside  the  furnace  to  minimize  external  oxidation  of  the  tube 
at  reaction  temperatures.  Both  graphite  boats  and  platinum  lined  alumina 
combustion  boats  have  been  used  satisfactorily  in  the  formation  of  ordered 
perovskite  fluorides. 

Preparation  of  pure  Cs2NaNdFg  has  never  been  satisfactory.  Early  in 
the  program,  before  the  argon  filled  glove  box  and  vacuum  line  were  available 
for  scimple  handling,  it  was  unclear  whether  pure  Cs2NaNdFg  did  not  form 
because  of  thermodynamic  considerations  or  whether  the  hygroscopic  samples 
deteriorated  so  quickly  that  good  x-ray  or  spectroscopic  data  could  not  be 
obtained.  Subsequent  work  with  adequate  facilities  has  shown  the  former  to 
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be  the  case,  in  that  anhydrous  samples  continued  to  exhibit  two  Nd 
in  fluorescence  spectra  and  the  x-ray  powder  diffraction  patterns  always 
revealed  small  amounts  of  impurities  were  present. 

Aleonard  and  Pouzet’  have  indicated  that  the  larger  rare  earth  ions  do 

not  form  cubic  perovskites  of  the  CS2NaMFg  composition.  Accordingly,  a 

preparation  of  Cs2NaYFg  was  made  in  the  same  manner  as  the  Cs2NaNdFg  samples 

previously  indicated.  The  product  was  nonhygroscopic  and  x-ray  data  agreed 

with  those  listed  for  Cs^NaYF,  in  Ref.  9.  This  result  is  excellent  evidence 

^ ° 3+  3+ 

that  the  difference  in  size  of  Nd  (0.107  nm)  and  Y (0.095  nm)  is  an 
important  factor  in  the  formation  of  stable  ordered  perovskite  fluorides. 

Subsequent  preparations  of  a range  of  compositions  in  the  Cs2NaNd^Y^_^Fg 
system  have  shown  that  the  hygroscopic  nature  of  the  product  varies  with 
the  Nd  content.  Powdered  samples  with  x = 0 . 5 or  less  present  no  particular 
handling  problems  in  air,  but  extended  exposure  of  powdered  samples  with  x 
greater  than  0.25  results  in  some  deterioration.  One  or  more  samples  have 
been  prepared  at  each  of  the  compositions  with  x = 1,  0.75,  0.50,  0.25, 

0.10,  0.05,  0.01  and  0.001.  Each,  with  the  exception  of  the  pure  Nd  material 
discussed  earlier,  has  the  ordered  perovskite  structure  as  determined  by 
x-ray  powder  diffraction  data.  These  samples  have  been  examined  in  some 
detail  spectroscopically  and  the  experiments  are  reported  in  the  Spectroscopic 
Results  section. 

800  °C 

REACTION:  2CsF  + NaF  + NdFs  CsoJ'JaNdFe 

H r , A r ^ 

INSULATION  FURNACE  CORE 


WITH  REACTANTS 


Figure  3.  Preparation  of  Perovskite  Fluorides 
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NaNd  Y,  F.  and  NaEu  Y,  F. 

X 1-x  4 X 1-x  4 

Polycrystalline  samples  of  these  materials  were  made  by  the  same 
procedure  used  for  Cs2NaNd^Y^_^Fg , except  that  somewhat  higher  temperatures 
(900°  to  950°C)  were  required.  The  NaNd  Y,  F.  and  NaEu  Y.  _ F.  series  were 

X 4 X X X 4 

prepared  at  compositions  with  x = 1,  0.25,  and  0.01  for  spectroscopic 
studies  in  order  to  trace  the  identity  of  impurities  in  Cs2NaNd^Yj^_^Fg 
materials . 

Mixed  Halides 

The  mixed  halide  materials,  Cs~NaNd  Y.  F,  Cl  and  CsTNaNd  Y.  Br,  Cl  , 

2 X 1-x  6-y  y 2 x 1-x  6-y  y 

were  specific  attempts  to  reduce  the  site  symmetry  of  the  Nd3+  ion.  The 
preparation  of  fluoride-chloride  samples  was  tried  under  an  argon  flow  in 
the  manner  described  for  the  Cs2NaNd^Yj^_^Fg  system  (except  no  HF  flow  was 
used) . Later  when  high-density  graphite  containers  with  tight-fitting  lids 
became  available,  the  attempted  preparation  from  component  fluorides  and 
chlorides  was  repeated.  A third  method  simply  involved  heating  the  com- 
ponents for  reaction  in  an  evacuated  quartz  ampule.  The  result  of  this 
last  method  was  extensive  quartz  attack.  The  other  saimples  were  shown  to 
be  impure  by  x-ray  data.  Fluorescence  spectra  were  characteristic  of 
CS2NaNdClg  with  no  reduction  of  lifetime.  Despite  a number  of  attempts, 
including  slowly  cooling  (3°/hr)  a sample  in  a closed  high  purity  graphite 
container,  no  preparation  of  a mixed  fluoride-chloride  ordered  perovskite 
was  successful. 

The  preparations  of  mixed  bromide-chlorides  were  made  in  evacuated 
quartz  ampules  by  heating  the  components  at  850°  for  several  hours  and  then 
cooling  the  samples  at  perhaps  100°/hour.  The  x-ray  data  were  correct  for 
ordered  perovskites  and  the  lattice  parameters  were  suitably  intermediate 
between  those  of  Cs2NaNdClg  and  that  expected  for  CS2NaNdBrg.  However, 
as  indicated  in  the  Spectroscopy  Results  section,  the  emission  spectra 
consisted  of  rather  broad  bands  and  were  essentially  the  same  as  spectra 
for  the  pure  chloride.  It  is  likely  that  this  resulted  in  a range  of 
bromide-chloride  Nd^"*"  environments. 

Crystal  Growth  Experiments 

A number  of  attempts  were  made  to  grow  crystals  of  Cs2NaNdFg  by  the 
gradient  freeze  method  in  a zoned  furnace  as  described  for  Cs2NaNdClg 
growth.  In  this  case,  a special  graphite  of  high  purity  and  density 
(Poco  DFP-3-2*)  was  used  to  construct  the  closed  containers  for  both 

*Union  Poco  Graphite,  Inc.,  Decater,  Texas. 
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horizontal  and  vertical  use  shown  in  Figure  4.  The  graphite  containers 
were  cleaned  prior  to  use  in  a stream  of  argon  and  HCl  at  900°C  for  several 
hours,  and  then  transferred  to  an  argon-filled  glove  box  for  sample  loading. 
The  closed  containers  were  effective  in  limiting  CsF  loss  even  when  subjected 
to  temperatures  above  800°C  for  several  days.  Losses  were  typically  a few 
tenths  of  a gram  from  a sample  charge  of  7 to  8 grams.  Temperature  zones 
of  from  40°  to  about  150°C  over  the  4 in.  container  length  were  tried  at  a 
cooling  rate  of  about  3°C/hr.  Perhaps  the  most  crystalline  samples  were 
obtained  with  the  higher  gradients,  but  nothing  approaching  single  crystal- 
line material  was  obtained.  The  x-ray  powder  data  always  indicated  small 
amounts  of  impurities  (nonperovskite  phase)  to  be  present.  This  was  con- 
firmed in  the  spectroscopic  data  because  the  hexagonal  form  of  NaNdF^  was 
always  present.  It  seems  likely  that  Cs2NaNd^Yj^_^Fg  melts  incongruentiy . 

In  retrospect,  perhaps  single  crystals  might  be  obtained  using  CsF  or  CsF-NaF 
eutectic  as  a solvent  for  normal  flux  growth  procedures. 


CLOSED  GRAPHITE  CONTAINERS 


~ r= 

TOP  SIDE 

HORIZONTAL  BOAT 


B. 

VERTICAL 

CYLINDER 


CRUCIBLE 

TYPE 


Figure  4.  Containers  Used  for  Fluoride  Crystal  Growth 
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2.3  SPECTROSCOPIC  RESULTS 


CS2NaNdClg 

Emission  and  absorption  spectra  of  this  material  have  been  reported  by 
Tofield  and  Weber®  who  interpreted  them  as  magnetic  dipole  and  forced  electric 
(vibronic)  transitions.  Although  there  is  little  doubt  that  the  environment 

represents  a true  0,  symmetry  ( see  for  instance,  the  work  on  the  corresponding 

+ 3 ” +3 

Eu  compound  by  Serra  and  Thompson  (Ref.  10)1  , the  Nd  spectrum  is  far  more 

complicated  and  the  vibronic  contributions  make  the  assignments  somewhat 

tentative.  Using  an  Xe  arc  as  excitation  the  spectrum  appears  to  be  identical 

to  that  reported  in  Ref.  5 (Figure  5 top  trace*).  For  comparison  purposes, 

we  also  give  an  emission  spectrum  of  the  same  sample  excited  by  the  dye 

laser.  The  differences  are  indeed  minor  except  for  the  change  of  intensity 

in  the  shoulders  of  the  dominant  line  in  the  resonance  transition  group 

(850  nm  to  950  nm) . 

C32NaNdFg 

As  noted  in  Section  2.2,  this  material  is  not  stable,  is  highly  hygro- 
scopic and  over  a period  of  time  in  air  loses  the  appearance  of  crystallinity. 
Figure  6 shows  two  emission  spectra  taken  with  an  Xe  arc  excitation.  The 
tcp  trace  has  been  obtained  on  a fresh  sample  and  has  a rather  poorly  re- 
solved structure.  The  decay  time  of  the  luminescence  is,  however,  1.28  ms 
characteristic  of  a site  with  inversion  symmetry.  The  lower  trace  was  taken 
after  a few  weeks  exposure  of  the  sample  to  laboratory  air.  We  notice  a 
radical  change  in  the  spectrum,  accompanied  by  a decrease  in  intensity.  The 
lifetime  of  this  sample  has  decreased  to  2.5  us. 


Cs^NaY,  Nd  F, 
2 1-x  X 6 


Depending  on  the  value  of  x,  emission  spectra  and  decay  times  can  be 
divided  into  fairly  distinct  categories  which  we  shall  designate  as  follows: 


Type  A:  0.1  <x  <0.75 

Type  B:  0.01  <x  <0.1 

Type  C:  x <0.01 

For  X <0.75  the  compounds  are  considerably  more  stable  than  for  x = 1 and 
can  be  left  for  weeks  without  appreciable  change.  Eventually  Type  A com- 
pounds do  undergo  some  change  in  the  same  direction  as  the  undiluted 


♦Excitation  conditions  and  observed  decay  times  (t.)  are  given  on  the  right 
side  of  the  spectra. 
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Figure  5.  Room  Temperature  Spectra  of  Pure  Perovskite  Chloride  Material  Under 
Broadband  and  Monochromatic  Excitation 


compound,  but  the  process  is  very  much  slower.  An  immediate  characteristic 
of  the  Type  A compounds  is  the  appearance  of  two  different  lifetimes  de- 
pending upon  the  excitation  wavelength.  In  Figure  7,  we  see  two  oscilloscope 
decay  traces  taken  under  identical  conditions  except  for  the  wavelength 
setting  of  the  dye  laser.  The  top  trace  shows  a very  rapid  decay  from  a 
high  peak  level  whereas  the  bottom  trace  shows  only  a low  intensity  rapid 
component  followed  by  a long  lifetime  component.  The  two  lifetimes  are 
2.4  us  and  2.2  ms,  respectively. 


8M 


MO 


«obo 


1060 


Figure  6.  Room  Temperature  Emission  on  Spectra  (Broadband  Excitation).  The 
Top  Trace  is  for  a Fresh  Sample.  The  bottom  trace  represents  the 
same  sample  after  a few  weeks. 


Cs2NaYo  sNdp  sFg 


EXCITATION 


577  nm 


597  nm 


Figure  7.  Decay  Current  for  Two  Different  Excitation  Wavelengths 
Showing  the  Existence  of  Two  Lifetimes 
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Representative  emission  of  Type  A compounds  are  shown  on  Figure  8. 

The  top  trace  is  obtained  with  Xe  arc  excitation  and  shows  a multitude  of 
lines.  A greatly  simplified  spectrum  can  be  obtained  by  tuning  the  dye 
laser  to  a wavelength  which  maximizes  the  fast  or  the  slow  component  of  the 
decay.  A further  simplification  resulting  in  the  two  lower  traces  is 
obtained  with  the  boxcar  whose  gate  is  adjusted  to  time  resolve  either  the 
fast  or  the  slow  component  of  the  delay.  By  judicious  choice  of  excitation 
wavelengths  and  gates,  one  can  almost  completely  eliminate  any  overlap  of 
lines  in  the  two  spectra.  Overall,  however,  the  Xenon  excited  spectrum  is 
dominated  by  the  long-lifetime  species  especially  in  the  resonance  region 
(note  the  similarity  of  the  47  and  28  trace  in  the  range  of  850  nm  to  950  nm) . 
Upon  closer  examination  the  Xe  excited  trace  represents  a summation  of  the 
two  time-resolved  traces,  not  in  the  sense  of  absolute  intensity,  but  the 
presence  of  the  respective  lines  which  are  tabulated  in  Tables  2 and  3. 

Upon  comparing  Figures  6 and  8,  a striking  observation  can  be  made. 

The  short  lifetime  spectrum  of  the  diluted  samples  (middle  trace  of  Figure  8) 
appears  to  be  identical  to  the  spectrum  of  a deteriorated  undiluted  sample 
(bottom  trace  of  Figure  6) . 


TABLE  2 

SHORT  LIFETIME  LINES  (nm) * 


4 4 

p -►  T 

*3/2  *9/2 

4 4 

P “►  T 

*3/2  *11/2 

896.0  (L?) 

1065.0  VW 

893.6  S 

1062.4  W 

890.4  S 

1060.4  W 

888.0  (L?) 

1058.8  S 

880.0  S 

1057.0  S 

877.0  S 

1055.0  S 

874.5  W 

1053.0  VW 

870.4  W 

1050.6  W 

864.6  S 

1048.4  S 

863.2  W 

1047.4  S 

861.6  S 

1043.8  W 

859.0  S 

1039.8  W 

10  (12?)  Lines 

12  Lines 

*Obtained  from  spectra  taken  at  liquid 
nitrogen  temperature 

S - Strong 

W - Weak 

VW  - Very  Weak 

(L?)  - Possibly  belonging  to  long 
lifetime  spectrum 
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Figure  8.  Emission  Spectra  of  50%  Nd  Perovskite  Fluoride.  The  top  trace  was  obtained  with 
broadband  excitation.  Lower  traces  were  time  resolved  (TR)  and  employed 


TABLE  3 


LONG  LIFETIME  LINES  (nm)* 

(CS2NaYo_75Ndo^25^6’ 


4 4 

3/2  ^9/2 

4 4 

3/2  -^11/2 

950.0  S 

1058.0  S 

930.5  S 

1050.0  W 

928.0  S 

1046.0  W 

918.5  W 

1039.0  W 

909.0  W 

1032.0  S 

896.0  S 

1023.0  VW 

888.0  S 

1020.0  VW 

879.0  VS 

876.0  W 

872.0  W 

867.0  W 

856.5  W 

12  Lines 

7 Lines 

♦Obtained  from  liquid  nitrogen  temperature 


spectra,  except  some  lines  in 


transitions  which  were  too  weak  an^  had  to 
be  obtained  at  room  temperature 


11/2 


VS  - Very  Strong 
S - Strong 
W - Weak 
VW  - Very  Weak 


When  the  concentration  of  Nd  drops  below  10  atom  %,  a significant  change 
in  the  spectra  is  manifest.  This  is  illustrated  in  Figure  9,  The  Xenon 
trace  differs  now  very  little  from  the  monochromatically  excited  traces. 

The  spectrum  begins  to  be  dominated  by  the  short  lifetime  species. 

By  scanning  the  dye-laser  excitation  and  observing  the  shape  of  the 
decay  curves,  it  is  easy  to  find  the  excitation  wavelengths  at  which  the 
fast  and  slow  decay  portions  are  respectively  maximized.  The  intensity 
(specifically  of  the  slow  component)  is,  however,  too  small  to  obtain  time- 
resolved  spectra  with  the  boxcar.  In  each  case,  therefore,  the  overall 
shape  of  the  spectrum  is  dominated  by  the  fast  component.  It  should  be 
noted  that  for  this  concentration  the  short  lifetime  has  lengthened  12  ys 
while  the  long  lifetime  has  shortened  to  0.9  ms.  This  represents  a truly 
intermediate  case.  When  the  concentration  of  Nd  is  dropped  to  1%  or  below 
(Type  C) , any  vestiges  of  two  centers  disappear.  This  is  illustrated  on 
Figure  10  where  broadband  and  monochromatic  excitation  produce  essentially 
the  same  spectrum.  Moreover,  only  one  lifetime  is  observed  (90  ys  for  this 
case)  and  there  is  no  detectable  change  of  spectrum  with  varying  excitation 
wavelength. 
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Figure  9.  Emission  Spectra  of  5%  Nd  Perovskite  Fluoride  Under  Broadband  Excitation  (Top  Trace) 

and  Two  Monochromatic  Excitations  Chosen  to  Maximize  Short  and  Long  Lifetime  Components 


Figure  10.  Emission 


At  a still  lower  concentration  (0.1%  Nd)  the  spectra  remain  essentially 
unchanged  but  the  lifetime  reaches  300  us.  The  diluted  samples  are  thus 
characterized  by  a single  center. 

Mixed  Hexahalides  Cs.,NaNdCl,  H 

£ °~y  y 

Many  attempts  were  made  to  remove  the  center  of  inversion  at  the  Nd 
site  by  substitution  of  mixed  halides  in  the  coordination  sphere  of  Nd. 

These  results  were  not  successful  as  shovm  on  Figure  11  which  cor’pares 
spectra  of  the  pure  chloride  with  that  of  mixed  chloride-bromide.  The 
effect  of  substitution  is  primarily  a broadening  of  the  spectrum  without 
drastic  effects  on  intensity  distribution  or  lifetime.  Mixed  chloride- 
fluoride  materials  did  not  accomplish  this  goal  either.  Substitution  of 
even  one  fluorine  atom  resulted  in  the  immediate  appearance  of  a very  short 
lifetime  site  shown  on  Figure  12,  whereas  the  bulk  of  the  emission  consisted 
of  essentially  the  unchanged  hexachloride  spectrum  (Figure  12,  top  trace) . 

NaNd  Y.  F. 

X 1-x  4 

As  explained  in  Section  2.2,  this  material  is  not  expected  to  yield 
an  efficient  stoichiometric  compound.  Our  interest  in  it  results  from  a 
remarkable  resemblance  of  its  emission  spectrum  with  the  short  lifetime 
emission  spectrum  of  the  Cs2NaNd^yj^_^Fg  compounds.  A direct  comparison  of 
the  two  spectra  is  shown  on  Figure  13 , the  top  trace  being  a reproduction 
of  the  middle  trace  (#29)  of  Figure  8. 

At  greater  dilution  the  tetraf luoride  spectrum  sharpens  and  again 
becomes  identical  to  the  dilute  hexafluoride  spectrum  shown  on  Figure  10. 
These  results  indicate  that  the  origin  of  the  short  lifetime  species  in  the 
hexafluoride  is  an  inclusion  of  the  tetraf luoride  phase.  This  latter  phase 
becomes  dominant  at  low  concentration.  Moreover,  since  it  is  apparently 
more  stable,  it  is  the  phase  into  which  the  hexafluoride  is  converted  by 
exposure  to  moisture  (see  Figure  6).  However,  even  Cs2NaNd^Yj^_^Fg , which 
were  carefully  prepared  and  handled  under  inert  conditions,  exhibited  the 
presence  of  the  short  lifetime  species  which  is  characteristic  of  NaNdF^. 

Eu^^  Compounds 

To  confirm  the  presence  of  a tetraf luoride  phase  in  the  fluoride 
elpasolites,  a series  of  analogous  compounds  containing  Eu^^  instead  of  Nd^^ 
were  prepared.  Since  the  level  structure  of  the  Eu^^  ion  is  considerably 
simpler  than  that  of  Nd^^,  the  spectra  are  simpler  to  interpret.  Samples 
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Figure  11.  Comparison  of  the  Emission  Spectra  of  Cs2NaNdCl6  (Bottom  Trace)  with  those 
of  Mixed  Chloride-Bromides.  Note  the  long  lifetimes  but  broadened  nature  i 
the  latter. 


Figure  13.  Direct  Comparison  of  the  Short  Lifetime  Spectra  Found  in  Hexafluoride  (Top  Trace) 
with  that  of  a Tetraf luoride  (Bottom  Trace) 


7 5 

were  pumped  with  the  dye  laser  tuned  to  the  F.-  D,  transition  occurring 

5 5 

in  the  vicinity  of  525  nm  and  and  emissions  occurring  from  the  Dq  or 
levels  to  the  lower  members  of  the  manifold  were  recorded. 

Figure  14  gives  a comparison  of  three  relevant  spectra.  The  two  upper 
traces  obtained  on  the  same  sample,  but  differing  only  in  the  excitation 
wavelength,  revealed  the  existence  of  two  different  sites  in  the  elpasolite. 

The  top  spectrum  is  an  excellent  example  of  an  octahedral  site  where  the 

5 7 5 7 

dominant  transitions  Dq  - F^  and  Dq  - F^  are  of  magnetic  character. 
Changing  the  excitation  from  525.5  nm  to  524.8  nm  results  in  a complete 
change  o^  the  spectrum  (middle  trace)  indicating  a lowered  symmetry.  The 
bottom  trace  is  obtained  from  a tetraf luoride  (NdEuQ  25^0  75^4^  containing 
the  same  molar  fraction  of  Eu  as  the  hexafluoride.  It  is  seen  to  be  very 
nearly  a replica  of  the  middle  trace.  This  suggests  that  here  also  the 
"distorted"  site  in  the  hexafluoride  corresponds  to  a site  characteristic 
of  the  tetraf luoride. 

An  interesting  feature  of  the  spectra  corresponding  to  the  tetraf luoride 

site  (Figure  14,  E8 , E19)  is  the  amost  complete  absence  of  lines  originating 

from  the  ^D^  state  when  time  integrated  detection  is  used.  This  indicates 
5 5 

a fast  Dj^  - Dq  relaxation.  When  the  concentration  of  Eu  is  reduced  to 
1%,  the  lines  originating  from  the  ^D^^  state  begin  to  appear.  This  is 
shown  on  Figure  15,  which  is  analogous  to  Figure  14  except  for  greater 

dilution  of  the  Eu  ions.  The  ^D,  - ^Dq  relaxation  is  thus  strongly  con- 

5 5 

centration  dependent.  This  suggests  that  it  may  occur  via  a D,  - D- : 

7 7 10 
Fq  - Fj  cross  relaxation,  which  would  obviously  depend  upon  the  average 

distance  between  the  ions.  It  is,  we  believe,  quite  significant  (see 

Discussion  section)  that  the  appearance  or  disappearance  of  lines  originating 

from  ^Dj^  levels  is  quite  parallel  in  the  pure  tetraf luoride  and  tetraf luoride- 

li)ce  site  in  hexafluoride. 

5 5 

The  dynamics  of  Dj^  - Dq  relaxation  were  studied  in  greater  detail  by 

using  methods  of  time-resolved  spectroscopy.  As  an  example  we  show  in 

Figure  16  time-resolved  traces  of  the  emission  ta)cen  at  two  different  times 

5 

during  the  decay.  At  early  times  (delay  0.5  ms, trace  E22)  the  D^^  emissions 
are  clearly  visible  whereas  at  longer  delays  (trace  E23)  they  are  absent. 

The  time  integrated  spectrum  of  the  same  material  (Figure  13,  trace  E19) 

shows  only  a very  wea)<  trace  of  the  ^D.  emission.  The  relaxation  rate  for 
5 5 ^ 

the  Dj^  - Dq  process  is  thus  relatively  slow,  of  the  order  of  a millisecond, 
at  a 25  atom  % concentration  of  Eu^'*’. 
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Figure  15.  Emission  Spectra  at  85'^K  of  the  Two  Eu'^  Sites  in  Cs2NaEu^  01^.  99^6  (Top  and  Middle 
Trace)  Compared  with  the  Site  in  Tetraf luoride  NaEu  (Bottom  Trace) . 
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There  is  one  further  remark  which  has  to  be  made  concerning  the  tetra- 
fluoride  spectra.  Invariably  in  the  region  we  observe  four  lines 

while  only  three  should  be  present  if  maximum  splitting  occurs.  The  same 
is  true  in  trace  E9  of  Figure  14  which  corresponds  to  the  "impurity"  site 
in  the  hexafluoride.  This  is  presumably  caused  by  the  presence  of  two 
different  sites  in  the  tetraf luoride  structure. 

Both  NaNd  y,_  F.  and  NaEu  Y,  F.  are  isostructural  with  the  hexagonal 
form  of  NaCeF^  which  has  the  6-Na2ThFg  structure.'*  In  this  structure 
2/3  of  the  Ce  atoms  occupy  nine  coordinate  sites  and  1/3  of  them  are  found 
in  six  coordination,  which  is  consistent  with  our  observation  that  Eu^^ 
occupies  two  sites  in  the  tetraf luoride.  We  have  not  been  able  to  resolve 

these  two  sites  at  liquid  nitrogen  temperature  by  fine  tuning  the  dye  laser. 

7 5 

Presumably  the  splitting  of  the  Fq  - transitions  (corresponding  to  the 

524.8  nm  excitation)  in  these  two  sites  are  separated  by  less  than  the 
minimum  bandwidth  of  our  dye  laser  source. 

Thus,  to  sum  up,  the  results  obtained  with  europium  compounds  support 
our  conclusion  that  the  short  lifetime  species  always  present  in  Nd  hexa- 
fluorides is  due  to  a contamination  by  the  tetraf luoride  phase. 

Lifetimes 

The  decay  times  of  the  Cs2NaNdjjY^_^Fg  compounds  show  the  existence 
of  two  vastly  different  values.  The  long  lifetime  is  characteristic  of 
the  hexafluoride  (cubic)  spectrum  while  the  short  lifetime  corresponds  to 
the  tetraf luoride  impurity.  Both  of  these  change  with  the  Nd  concentration 
and  have  been  plotted  on  the  same  (Figure  17)  graph.  Both  lifetimes  show 
an  increase  with  dilution,  the  long  lifetime  reaching  a remar)cably  large 
value  of  5 ms  at  x = 0.1.  The  short  lifetime  increases  monotonically  from 
about  3 US  at  X = 1 to  300  ms  (point  not  included  on  the  graph)  at  x = 0.01. 
Below  X = 0.1  (roughly  the  demarcation  between  A and  B type  spectra)  the 
intensity  of  the  long  lifetime  emission  becomes  too  weak  to  obtain  reliable 
measurements  of  the  decay. 

We  should  point  out  here  that  the  x-ray  powder  data  show  the  dominant 
crystalline  phase  present  to  be  the  ordered  perovskite;  this  is  especially 
so  for  samples  with  a low  Nd  content.  In  fact,  in  many  samples,  it  is  dif- 
ficult to  ascertain  the  presence  of  the  tetraf luoride  phase  from  x-ray  data. 
It  seems  reasonable  to  suppose  the  short  lifetime  Nd  emission  is  more  intense 
than  that  of  the  long  lifetime  site.  The  reduced  long  lifetime  emission  at 
concentrations  below  x = 0.1  probably  reflects,  in  part,  the  reduced  Nd 
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short  lifetime  (la) 


3.  DISCUSSION 


3 . 1 MATERIALS 

Very  Little  needs  to  be  added  to  the  description  already  given  of 
Cs2NaNdClg  and  related  materials.  Both  the  x-ray  and  spectroscopic  data 
agreed  that  the  Nd  ion  is  present  in  a single  site  of  Ojj  symmetry.®'^'*” 

One  might  only  add  that  sections  of  S2unple  rods  were  difficult  to  polish 
because  of  the  softness  and  moisture  sensitivity  of  the  chloride  material. 
This  would  seem  to  severely  limit  the  utility  of  perovskite  chlorides 
as  potential  laser  materials. 

Most  of  this  work  encompassed  the  investigation  of  the  Cs2NaNd^Yj^_^Fg 
system.  The  obvious  conclusion  is  that  there  is  a strong  tendency  toward 
a partial  decomposition  to  CsF  and  NaNdF^.  Several  groups  had  reported 
the  preparation  of  Ca2NaYFg  by  reaction  of  pressed  pellets  of  mixed  com- 
ponents at  500°  to  600°C.’''^  No  mention  of  decomposition  was  given  at 
higher  temperatures  and  it  is  presumed  that  the  principal  reasons  for  the 
use  of  such  low  temperatures  were  to  avoid  container  attack  and  to  minimize 
loss  of  the  volatile  CsF. 

From  our  work,  in  which  a number  of  samples  were  annealed  overnight 
below  600°C,  it  seems  unlikely  that  really  pure  phase  material  results  even 
from  low  temperature  preparations.  It  would  be  helpful  to  obtain  detailed 
phase  diagram  data  for  this  system,  but  this  was  beyond  the  scope  of  this 
project.  Since  CsF  is  both  volatile  and  insensitive,  it  would  be  necessary 
to  design  closed  containers  suitable  for  thermal  analyses.  Of  course, 
special  care  in  handling  the  materials  would  be  required  whether  the  system 
were  investigated  as  the  ternary  (CsF,  NaF,  NdF^)  or  the  pseudobinary 
(CsF,  NaNdF^). 

As  has  been  pointed  out,  attempts  to  distort  the  Nd  site  in  ordered 
perovskites  were  unsuccessful.  These  primarily  involved  anion  substitution 
such  as  the  partial  substitution  of  bromide  for  chloride.  The  chief  result 
was  a broadening  of  the  emission  spectra,  presumably  from  the  presence  of 
multiple  Nd  sites  from  random  substitution.  The  attempted  substitution  of 
fluoride  for  chloride  failed  to  form  a genuine  mixed  anion  perovskite.  It 
would  seem  that  for  a single  site  to  result  that  ordering  of  the  anions 
would  be  necessary.  This  would  probably  imply  the  necessity  of  a significant 
difference  in  anion  size  and/or  a difference  in  anion  charge.  Perhaps  oxy- 
fluorides  or  oxychlorides  might  offer  a better  prospect  for  a reduction  of 
Nd  site  symmetry. 
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Variation  of  M in  the  formula  Cs^MNdF,  was  also  tried  with  M represented 

J o 

by  Li  and  K,  and  mixtures  of  Mg  and  Na.  The  potassium  substitution  seemed 
superficially  to  give  a cleaner  x-ray  powder  pattern,  but  the  superlattice 
lines  were  very  much  reduced  in  intensity  and  a short  lifetime  site  present 
in  the  fluorescence  spectra  was  assigned  to  that  in  KNdF^.  The  Li  sub- 
stitution failed  to  result  in  a perovskite  material  and  the  attempted  dis- 
tortions with  Md  substitutions  also  were  not  useful. 

3 . 2 SPECTROSCOPY 

The  spectra  of  all  the  materials  studied  by  us  depart  significantly 
from  the  simple  picture  of  a Nd'^  ion  placed  in  a site  of  octahedral  sym- 
metry. The  expected  transitions  allowed  by  symmetry  selection  rules  are 
shown  on  Figure  18 . ' ’ As  shown,  the  magnetic  dipole  transitions  would 

result  in  three  lines  in  the  resonance  region  (^F-/-  -*•  ^I-,_)and  four  lines 

l4  4 \ ^ 3/^' 

in  the  laser  transition  region  ^ ^3/2  ^11/2/  * departure  from  cubic 

symmetry  removes  all  the  J + 1/2  degeneracy  (except  Kramer's  degeneracy) 

4 4 4 

from  all  the  levels  resulting  in  splitting  of  the  ^3/2'  ^11/2  ^9/2 

levels  into  respectively  two,  six  and  five  Stark  components.  This  should 
cause  the  appearance  of  12  lines  in  the  laser  transition  region  and  10  lines 

in  the  resonance  region  as  shown  on  Figure  19  for  the  case  of  D..  symmetry 

3+ 

(representing  a tetragonal  distortion  of  the  octahedral  Nd  site) . The 
removal  of  all  symmetry  forbidden  transitions  in  Nd^"*^  upon  lowering  the 
symmetry  makes  the  interpretation  of  the  spectra  of  that  ion  much  more  dif- 
ficult than,  say,  Eu^^.  Serra  and  Thompson  report,*”  for  instance,  a strict 
adherence  to  the  0^  symmetry  selection  rules  in  the  compounds  Cs2NaEuClg 
and  Cs2Na(Eu,  Y)Clg.  The  case  of  Nd  compounds  is  evidently  cjuite  different 
since  invariably  more  than  a total  of  7 lines  (characteristic  of  Oj^)  and 
often  more  than  22  lines  (symmetry  less  than  Oj^)  are  observed.  This  leaves 
open  the  possibilities  of  vibronic  transitions  and  multiple  sites  or  com- 
binations of  both. 

The  results  presented  in  Section  2 indicate  that  in  the  majority  of 
fluoride  compounds  at  least  two  sites  are  present.  Referring  to  Table  2, 
we  see  that  the  short  lifetime  site  can  be  interpreted  as  due  to  a dis- 
tortion (lower  symmetry  than  Oj^)  without  the  need  to  invoke  vibronic  transi- 
tions. If  the  896.0  and  888.0  nm  lines  are  rejected,  the  total  number  of 
10  and  12  lines  is  consistent  with  the  2 (J  + 1/2)  of  possible  transitions. 
This  is  consistent  with  the  interpretation  of  the  short  lifetime  site  as 
due  to  the  tetraf luoride  phase. 
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4 3 


Figure  18.  Level  Splitting  and  Symmetry  Allowed  Transitions  in  Symmetry 

The  situation  is  more  complicated  for  the  long  lifetime  spectrum 
(Table  2) , since  in  this  case,  we  observe  more  than  the  maximum  possible 
number  of  lines  (12  instead  of  10)  in  the  resonance  region  and  7 instead  of 
the  maximum  of  12  in  the  laser  transition  region.  This  suggests  that  the 
long  lifetime  site  may  have  nearly  octahedral  symmetry  but  that  some  of  the 
lines  have  vibronic  origins.  This  is  the  conclusion  reached  by  Tofield  and 
Weber  in  their  work  on  Ca2NaNdClg.® 

The  interpretation  of  vibronic  spectra  is  usually  quite  a complicated 
matter  (especially  for  an  ion  like  Nd^^  showing  a multitude  of  lines) , 
requiring  emission  and  absorption  measurements  often  at  helium  temperatures. 
It  is  futile  to  base  interpretation  only  on  emission  spectra  taken  on  powder 
S4unple8. 

The  dependence  of  lifetime  upon  concentration  shown  in  Figure  17  in- 
dicates that  throughout  most  of  the  composition  range  the  luminescence  of 
both  the  long  and  short  lifetime  sites  is  strongly  concentration  quenched. 
The  laws  of  this  quenching  are,  of  course,  of  great  importance  if  better 
stoichiometric  materials  are  to  be  developed. 
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Figure  19.  Level  Splitting  and  Symmetry  Allowed  Transitions  in  Symmetry 

It  is  customary  to  express  the  quenching  process  by  the  following 
expression: 


^d 


lCc(x) 


where  is  the  observed  decay  of  luminescence,  the  decay  at  infinite 

dilution  and  k (x)  a rate  constant  for  the  quenching  process  which  is  a 
® 3+ 

function  of  x,  the  fraction  of  Nd  in  the  rare  earth  sites.  A normalized 
(dimensionless)  version  of  this  constant  is  given  by: 


T 


O 


k_(x) 


I 

It  has  been  noted  recently  by  Singh,  et  that  for  Nd  pentaphosphate 

the  quenching  constant  varies  linearly  with  the  Nd^^  concentration.  Chinn, 
Hong  and  Pierce' have  confirmed  this  observation  and  have  shown  that  the 
linear  relationship  is  obeyed  by  several  other  stoichiometric  materials 
studied  by  them.  The  same  appears  to  be  true  of  the  Cs2NaNd^Yj^_^Fg  system 
as  shown  on  Figure  20.  The  upper  curve  corresponds  to  the  short  lifetime 
species  (with  taken  as  300  us  corresponding  to  x = 0.01)  and  the  bottom 
curve  corresponds  to  the  long  lifetime  component.  In  this  case,  the  choice 
of  is  rendered  difficult  by  the  lack  of  data  below  x = 0.1  (see  Figure  17). 
We  have  chosen  to  be  9 ms  which  is  a rough  extrapolation  of  the  long 
lifetime  curve  to  x = 0 in  Figure  17.  This  behavior  is  in  marked  contrast 
to  doped  materials  such  as  Nd:YAG,  NdrLaF^;  Nd;La202S  in  which  the  quenching 
rate  varies  as  the  square  of  the  concentration.  It  is  interesting  to  note 
that  the  linear  quenching  rate  applies  also  to  a material  like  Cs2NaNdj^Y^_^Fg 
and  thus  to  situations  where  electric  dipole  transitions  (at  least  for  the 
long  lifetime  site)  may  be  forbidden. 


Figure  20.  Normalized  Rate  Constant  for  Concentration  Quenching  in 
CsjNaYj^.jjNdjfFg : •-Short  Lifetime  Luminescence;  x - Long 
Lifetime  Luminescence.  Data  from  Figure  16. 
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The  short  lifetime  center  occurring  in  Cs^NaNd  Y,  F,  has  been  ascribed 

to  the  presence  of  a luminescence  center  similar  to  that  occurring  in 

tetraf luoride.  There  is  at  present  no  definite  evidence  whether  the  tetra- 

fluoride  appears  as  a separate  phase  or  as  a localized  distortion  of  the 

lattice.  The  short  lifetime  center  is  much  more  strongly  quenched  than  the 

long  lifetime  center.  The  absolute  values  of  k for  the  short  lifetime 

cn 

centers  are  more  than  two  orders  of  magnitude  larger  than  those  reported 
for  other  stoichiometric  materials.  For  the  long  lifetime  centers,  the 
disparity  is  much  less  - only  about  a factor  of  3 or  4.  The  meaning  of  this 
is  not  entirely  clear  in  the  absence  of  better  structural  information  and 
better  understanding  how  quenching  occurs  in  the  first  place.  It  is  sig- 
nificant, however,  that  we  have  not  been  able  to  detect  any  energy  transfer 
between  the  two  systems  of  centers.  Each  can  be  excited  individually  by 
choosing  the  right  pxcitation  wavelength.  This  lack  of  communication  between 
centers  is  a strong  argument  for  the  existence  of  a separate  tetraf luoride 
phase  spatially  separated  from  the  elpasolite  phase.  This  conclusion  is 
also  supported  by  the  sensitivity  of  the  relaxation  in  Eu^'*’  com- 

pounds as  a function  of  concentration,  mentioned  in  Section  2.3.  If,  the 
tetraf luoride  sites  were  randomly  dispersed  in  the  hexafluoride  lattice 
one  would  expect  little  tendency  for  cooperative  phenomena  to  occur. 

The  work  reported  in  this  paper  indicates  that  the  design  of  a stoichio- 
metric Nd  material  based  on  the  elpasolite  structure  but  with  increased 
radiative  probability  (center  of  inversion  removed)  is  quite  a bit  more 
difficult  than  originally  anticipated. 

The  fluoride  materials  have  proved  to  be,  in  general,  less  stable  than 
the  chlorides  and  invariably  led  to  the  existence  of  two  centers.  Also, 
it  seems  unlikely  that  a distorted  perovskite  of  chloride  or  fluoride  would 
be  useful  as  a laser  material  if  Cs  is  involved  as  the  large  cation  because 
of  instability  in  moist  air.  Nevertheless,  this  work  has  led  to  some  con- 
clusions which  are  of  intrinsic  interest  and  are  also  of  value  for  future 
work  in  this  area. 
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SECTION  II 

POTASSIUM  LITHIUM  NEODYMIUM  FLUORIDES 
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1 . INTRODUCTION 


This  work  is  primarily  comprised  of  the  investigation  of  a particular 
Nd  containing  phase  in  the  KF-LiF-NdF^  system,  which  from  spectroscopic 
studies,  shows  promise  as  a mini-laser  material.  The  investigation  of 
various  compositions  in  the  KF-LiF-NdF^  ternary  system  was  an  outgrowth 
of  our  investigation  of  neodymium  ordered  perovskite  halides  as  potential 
mini-laser  candidate  materials.  As  discussed  in  the  report  on  these  latter 
materials,*  a principal  problem  was  the  removal  of  the  center  of  inversion 
at  the  octahedral  neodymium  site.  A number  of  attempts  were  made  at  this 
which  primarily  involved  anion  and  cation  substitutions.  A third  approach 
was  to  exceed  the  limits  of  cubic  perovskite  stability  as  dictated  by 
the  relative  size  of  the  ions  of  the  A2BB'Xg  composition.  Among  oxides 
of  the  perovskite  type,  the  tolerance  factor,  t,  represents  the  situation 
where  all  the  component  ions  are  in  contact  and  is  equal  to  one  when  the 
A and  X ions  are  the  same  size  in  the  general  ABX^  cubic  perovskite.  In 
the  following  expression. 


t = ^ ^ 

/2  (Rg  t R^) 

R^,  Rg  and  R^^  refer  to  the  ionic  radi’ - ‘ For  cubic  perovskite  oxides, 
t falls  between  0.8  and  0.9.  A number  of  oxides  are  known  in  which  the 
unit  cell  has  been  distorted  to  tetragonal,  hexagonal,  orthorhombic  or 
monoclinic  types.  This  does  not  necessarily  mean  that  the  nearest 
neighbor  B ion  symmetry  has  been  changed  very  much,  if  at  all. 

Perovskite  halides  of  less  than  cubic  symmetry  are  somewhat  rarer. 
Hexagonal  and  tetragonal  distortion  have  been  observed  for  fluorides  as 
t becomes  less  than  "v-C.B  or  greater  than  ''^0.9.  Monoclinic  structures  have 
been  reported  for  KMgF^  (t  = 0.94),  RbMgF^  (t  = 1.03)  and  KCaF^  (t  = 0.93). 
The  expression  for  the  tolerance  factor  becomes  less  useful  for  ordered 
perovskite  halides  containing  rare  earth  ions  because  of  the  larger  size 
differential  (e.g.,  for  CS2NaNdFg)  between  the  A and  X ions  and  the  com- 
plicating factor  that  Rg  represents  two  ions  of  different  size  and  charge. 

This  is  probably  less  true  for  a compound  like  Cs-NaNdCl,  in  which  Cs^  and 

— + 3+  ^ ” 

Cl  are  about  the  same  size  and  Na  and  Nd  are  not  too  different  in  size. 


*See  Section  I 
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The  calculated  t for  this  material  is  0.87  to  0.89,  depending  on  the  value 
of  Rg  chosen,  which  is  well  within  the  range  for  perovskites. 

Certain  ions  tend  toward  nonsymmetric  octahedral  coordination,  such 
as  Cu^^  in  CsCuCl^,  in  which  two  of  the  Cl  are  more  remote  than  the 
other  four. ^ Octahedral  coordination  is  not  very  common  for  trivalent 
rare  earth  ions  which  are  mostly  found  in  seven  to  nine  coordinate  sites. 
Neodymium  is  certainly  found  in  a regular  octahedral  site  in  Cs2NaNdClg . " > ® 
In  general,  similar  fluorides  might  be  expected  to  show  less  tendency 
toward  regular  octahedral  rare  earth  sites  because  the  fluoride  ion  is 
much  smaller  than  the  chloride  ion. 

Preparations  of  a number  of  A2BB'Xg  ordered  cubic  perovskites  were 
made  in  the  course  of  this  project  (see  Section  I) . None  showed  x-ray 
evidence  of  any  distortion  of  the  cubic  unit  cell.  These  were  mostly 
fluorides  with  A ions  represented  by  the  larger  group  I cations  (Cs,  Rb,  K) , 
the  smaller  B ions  by  Li  through  Rb,  and  B'  ions  were  various  trivalent 
rare  earths  (especially  Nd) . Among  fluorides,  the  rather  large  Nd  ion  was 
found  to  lead  to  a partial  decomposition  of  the  perovskite  phase  to  the 
tetrafluoride  (NaNdF^,  KNdF^)  and  an  alkali  fluoride.  The  composition 
K2LiNdQ  25'“"o  75^6  among  those  prepared  and  examined  spectroscopically. 

Interestingly,  although  the  composition  contained  several  phases,  one  of 
them  gave  spectroscopic  indications  of  a long  lifetime  component,  which 
in  turn  implied  a large  separation  of  the  Nd  ions  in  the  responsible  phase. 
This  report,  then,  describes  our  investigation  of  the  KF-LiF-NdF^  system 
and  the  search  for  a new  mini-laser  material. 
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2.  MATERIALS  PREPARATION 

Polycrystalline  samples  for  the  potassium  lithium  neodymium  fluoride 
system  were  routinely  prepared  by  heating  the  anhydrous  component  fluorides 
in  a flow  of  argon  and  HF  at  about  900°C  for  two  or  three  hours.  The 
samples  were  located  in  a small  platinum  foil  boat,  which  was  placed  in 
a larger  alumina  combustion  boat  for  support.  The  Monel  apparatus  and 
the  general  procedure  which  was  used  was  that  described  elsewhere  (see 
Section  I -2.2)  for  the  preparation  of  ordered  perovsJcite  samples.  The  reac- 
tants were  stored  in  an  argon-filled  glove  box,  although  only  KF  is  partic- 
ularly moisture  sensitive.  Starting  materials  (Kf,  LiF,  NdF^)  were  dried  under 
dynamic  vacuum  at  100°  to  150°C.  However,  no  moisture  was  collected  in  a 
cold  (-195°C)  trap  and  no  differences  were  found  in  reaction  products  made 
from  reactants  which  were  not  dried  prior  to  use.  The  KF  and  LiF  were 
five  9'2  pure  from  United  Mineral  and  NdF^  was  supplied  in  four  9's  purity 
by  Research  Chemical.  Weight  losses  in  a typical  preparation  were  generally 
less  than  O.lg  from  2.0g  samples  made  in  a few  hours  under  flow  conditions. 

In  colosed  graphite  boats,  7g  samples  were  exposed  to  temperatures  over 
800°C  for  several  days  with  weight  losses  of  O.lg  to  0.2g.  As  a 
precautionary  measure,  reaction  products  were  handled  under  argon  in  a 
glove  box. 

Samples  of  the  following  compositions  were  prepared  for  studies  in 
this  program:  K^LiNdFg,  K^LiNdF^,  K^LiNdFg,  KLiNdF^,  KNdF^  and  K^NdFg. 

These  were  examined  by  x-ray  powder  diffraction  techniques  as  well  as  the 
special  methods  of  Section  3 of  this  Report. 

As  more  became  known  of  the  phase  of  interest,  several  series  of 
diluted  samples  were  synthesized.  These  included  K_LiNd  Y,  F,., 

/ X i.~X  b 

K2LiNd^Gd2^_j^Fg , K2LiNd^La^_^Fg  and  K2LiNd^Cej^_j^Fg . Of  these,  the  most 
useful  was  the  K2LiNd^Cej^_^Fg  because  it  provided  a wider  isostructural 
composition  range. 

A number  of  attempts  have  been  made  to  grow  single  cyrstals  of  the 
Nd  phase  of  interest.  The  gradient  freeze  method  has  been  employed  using 
the  closed  horizontal  and  vertical  graphite  containers  sketched  in  Figure  1. 
These  were  constructed  from  a special  high  density,  high  purity  graphite 
(Poco  DFP-3-2*).  Since  the  melting  point  is  uncertain,  samples  of 
K.,LiNdF,  and  K,LiNdF,  have  been  cooled  slowly  (2°  to  5°C/hr)  from  about 

Zb  J b 

850°C  with  gradients  of  25°  to  100°C  along  the  4-in.  boat.  The  graphite 
containers  were  maintained  in  an  argon  environment  during  the  crystal 
growth  process.  In  early  experiments  severe  carbon  contamination  was 

♦Union  Poco  Graphite,  Inc.,  Decater,  Texas. 
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Figure  1.  Closed  Graphite  Containers  Used  for  Crystal  Growth 

obvious  even  though  precautions  were  taken  to  dry  the  samples  prior  to 
growth  under  dynamic  vacuum.  On  occasion,  the  starting  materials  were 
pretreated  by  heating  in  Pt  containers  up  to  800°C  in  an  HF-Ar  gas  flow. 
The  situation  was  finally  improved  markedly  by  pretreating  the  graphite 
containers  in  a flow  of  HCl  at  950°C,  then  transferring  them  under  Ar  to 
an  inert  atmosphere  glove  box  where  sample  loading  took  place.  Even  this 
was  not  entirely  satisfactory  with  respect  to  graphite  contamination. 

Later  when  a flux  growth  method  (KF*2LiF)  was  employed  and  the 
container  shape  was  not  considered  important,  the  graphite  boats  were 
lined  with  one  mil  platinum  foil  so  as  to  further  reduce  container 
attack.  This  sample  was  cooled  at  2J5°C/hr  from  about  800°  to  500°C. 
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3.  RESULTS 


The  equipment  and  methods  employed  in  the  spectroscopic  investigation 
of/ the  KF-LiF-NdF^  system  has  been  described  in  detail  (Section  I - 2.1) 
and  will  not  be  repeated  here.  The  emission  spectra  of  three  different 
fluorides  of  Nd  are  given  in  Figure  2.  All  of  these  were  taken  with  a 
broadband  Xenon  arc  excitation.  The  top  trace  is  that  of  the  tetraf luoride 
KNdF^,  a compound  in  which  nine  coordinate  Nd^^  ions  exist,  and  having 
an  orthorhombic  structure  with  the  space  group  Pnma.^  As  an  efficient 
material,  the  tetraf luoride  is  of  little  interest  since  its  lifetime 
(''-5  ps)  is  strongly  quenched  at  stoichiometric  composition.  The  spectrum 
associated  with  the  tetraf luoride  phase  is  designated  as  Type  I.  The 
middle  trace  (Type  II)  belongs  to  the  composition  K^LiNdF.^  and  is  of  con- 
siderably greater  importance  because  of  a much  longer  decay  time  (280  ys) 
indicating  significantly  reduced  quenching.  The  bottom  trace  (Type  III) 
belongs  to  the  composition  K^NdF^.  Although  the  emission  spectrum  has  an 
attractive  simplicity  and  apparently  favorable  branching  ratio  (ratio  of 
integrated  emission  in  the  1.05  nm  region  to  that  at  0.89  nm) , the  decay 
time  is  again  short  - of  the  order  of  a few  microseconds. 
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Figure  2.  The  Emission  Spectra  of  Three  Different  Nd  Fluorides 
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Clearly  the  fluoride  of  composition  K^LiNdF^  with  its  characteristic 
Type  II  spectrum  is  of  greatest  interest.  In  fact,  the  decay  time  of 
280  ys  is  the  longest  known  for  any  noncentrosymmetric  stoichiometric  Nd 
compound.  The  lack  of  sensitivity  to  the  exact  proportion  of  p>otassium  to 
fluorine  is  illustrated  in  Figure  3,  which  shows  two  identical  (Type  II) 
spectra  for  materials  of  nominal  composition  K^LiNdF^  and  K2LiNdFg. 


HAVELENQTH  (MICRONtl 


Figure  3.  The  Type  II  Emission  Spectra  Characteristic  of 
KjLiNdF.^  and  K2LiNdFg  Compositions 

Some  preparations  of  K2LiNdFg  show  a slight  distortion  of  the  Xenon 
excited  spectrum  particularly  noticeable  in  the  long  wavelength  tail 
(1.063  nm)  and  as  a short  wavelength  shoulder  (1.056  nm)  in  the  1.058  nm 
line.  This  is  illustrated  in  the  top  trace  of  Figure  4.  Upon  excitation 
with  a monochromatic  (dye  laser)  source  tuned  to  maximize  the  emission,  it 
becomes  immediately  obvious  that  the  broadband  excited  spectrum  is  a super- 
position of  a Type  II  spectrum  (with  excitation  peak  at  428.4  nm)  and  a 
Type  I spectrum  (with  an  excitation  peak  at  427.4  nm) . This  is  shown  on 
the  two  lower  traces  of  Figure  4.  The  overall  emission  is,  however, 
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KjLiNdFj 


Figure  4.  Emission  Spectra  of  K2LiNdFg  with  Xenon  Arc  and 
Selected  Dye  Laser  Excitation 


dominated  by  the  Type  II  spectrum,  the  Type  I ( tetraf luoride)  being  present 
only  as  a small  component.  The  dominance  of  the  Type  II  spectrum  is  even 
more  pronounced  in  the  Xenon  spectra  of  K^LiNdF.^  and  K2LiNdFg  samples 
shown  in  Figure  3. 

The  effect  of  dilution  (in  this  case  with  Ce^^) on  the  spectra  of 
K2LiNd^Cej^_^Fg  is  shown  on  Figure  5.  Here  all  the  traces  were  taken  with 
broadband  Xenon  excitation.  The  clearly  visible,  overall  effect  is  that, 
as  the  Nd  content  is  reduced,  the  spectrum  shows  a definite  transformation 
from  Type  II  to  Type  I.  At  1%  composition  (x  = 0.01),  the  emission 
spectrum  is  virtually  indistinguishable  from  that  of  the  tetraf luoride 
phase.  Thus  at  the  high  concentration  end  lx  1) , the  emission  is 
dominated  by  the  Type  II  spectrum,  the  tetraf luoride  component  (Type  I) 
being  relatively  very  weak.  Exactly  the  opposite  is  true  at  the  other 
end  of  the  composition  scale  (x  'v  0.01)  where  the  Type  II  spectrum  is  the 
much  weaker  component.  Figure  6 shows  the  Type  II  spectrum  still  present 
in  this  diluted  material.  Note  that  the  ratio  of  the  peak  intensities 
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of  the  long  and  short  wavelengths  components  has  decreased  considerably 
from  the  case  of  the  undiluted  compound  (Figure  3) . We  shall  use  this 
fact  in  Section  4,  Discussion,  to  present  a qualitative  picture  of  the 
composition  of  the  material. 


Lifetimes 

In  order  to  measure  decay  times,  the  dye  laser  was  tuned  to  the  main 
excitation  peaks  of  the  two  types  of  spectra  (427.4  nm  for  Type  I,  428.4  nm 
for  Type  II).  The  decay  was  followed  for  approximately  two  orders  of  mag- 
nitude with  the  boxcar  operating  in  the  logarithmic  output  mode. 

The  dependence  of  the  decay  times  on  composition,  for  the  two  types 
of  spectra,  are  given  in  Figure  7.  Both  curves  show  a decrease  of  life- 
time with  concentration,  the  major  portion  of  the  decrease  being  confined 
to  rather  low  concentration.  In  discussing  the  significance  of  these  curves, 
one  should  note  that  the  variable  x is  defined  only  by  the  nominal  com- 
position of  the  compound  K2LiNd^Cej^_^Fg . Since  there  is  no  guarantee  that 
the  ratio  of  this  compound  to  the  tetraf luoride  phase  KNdF^  is  constant 
throughout  the  range  of  x,  the  abscissa  of  each  of  the  experimental  points 
may  be  in  error.  Because  of  this  uncertainty,  we  prefer  to  represent  the 
lifetime  quenching  of  the  important  Type  II  emission  as  made  of  segments 
rather  than  a monotonic  function. 


Figure  7. 


The  Composition  Dependence  of  the  Fluorescence 
Lifetimes  of  Types  I and  II  Nd  Sites 
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For  the  Type  II  emission,  the  decay  is  strictly  exponential  over  the 
whole  range  of  x as  shown  in  Figure  8.  This  is  not  true  of  the  decay  of 
the  Type  I emission  which  shows  departure  from  exponential  behavior  as 
illustrated  on  Figure  9.  Nonexponential  decay  is  encountered  often  in 
luminescence  but  since  in  the  present  case.  Type  I emission  is  an  unwanted 
impurity  component,  we  have  not  investigated  this  aspect  in  any  detail. 
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Figure  8.  The  Exponential  Decay  of  the  Type  II 
Neodymium  Species  at  Different 
Temperatures 
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Figure  9.  The  Fluorescence  Decay  Curve  for  a Typical  Type  I 
Neodymium  Site.  Note  nonexponential  behavior. 

X-ray  powder  diffraction  patterns  of  the  various  samples,  as  well  as 
spectroscopic  data,  were  obtained  as  the  materials  became  available.  A 
comparison  of  the  x-ray  powder  data  for  KNdF^,  K^LiNdFg  and  K^LiNdF^  are 
shown  in  Figure  10.  It  is  obvious  that  lines  are  in  common  in  all  three 
samples.  Those  lines  corresponding  to  tetraf luoride  lines  became  pro- 
gressively weaker  as  the  "extra"  lines  became  stronger  in  going  to  the 
KjLiNdF^  composition.  Neither  K2lii2NdF^  or  K^LiNdFg  show  further  develop- 
ment of  intermediate  x-ray  lines,  and  if  anything,  produce  x-ray  patterns 
which  are  more  like  KNdF^ . Ordinarily  one  would  assume  the  patterns  of 
K2LiNdFg  and  K^LiNdF.^  to  represent  KNdF^  plus  one  or  more  additional 
phases  in  lesser  amounts.  The  materials,  as  routinely  prepared  in  a 
single  day,  are  not  particularly  crystalline  in  appearance  and  it  was 
not  possible  to  pick  out  individual  crystals  for  separate  examination. 
However,  the  spectroscopic  evidence  (Spectroscopic  Results  section)  clearly 
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Figure  10.  Comparison  of  the  X-Ray  Powder  Diffraction  Patterns  of  KNdF^, 
K2l<iN<3Fg  and  K^LiNdF^.  Relative  intensity  is  plotted  vs.  20. 

indicated  the  long  lifetime  phase  to  be  predominant.  This,  plus  the  fact 
that  not  very  much  difference  was  seen  in  either  x-ray  or  spectroscopic 
data  as  KF  was  added  to  the  compositions,  lead  us  to  believe  that  the  in- 
creased isolation  of  Nd^^  was  achieved  by  an  expansion  of  the  orthorhombic 
KNdF^  unit  cell®  and  that  the  extra  lines  were  possibly  superlattice 
lines.  Accordingly,  an  effort  was  made  to  index  these  patterns  by  using 
an  available  program  capable  of  indexing  cubic,  hexagonal,  tetragonal  or 
orthorhombic  structures.  This  effort  was  not  successful,  indicating  that 
if  the  unit  cell  were  expanded  that  this  might  well  be  monoclinic  rather 
than  orthorhombic.  Of  course,  lines  from  an  impurity  phase  could  have 
interfered  in  the  indexing  procedure. 
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The  crystal  growth  experiments  with  K2LiNdFg  and  K^LiNdF^  compositions 
have  been  somewhat  disappointing.  The  later  experiments  gave  samples  which 
are  considerably  more  crystalline  in  appearance  than  single  day  prepara- 
tions, but  obviously  are  not  single  crystals  as  a whole.  Some  material 
from  these  samples  has  been  ground  lightly  and  attempts  have  been  made 
to  pick  out  single  crystal  pieces  under  a microscope  at  low  magnification. 
Thus  far,  most  of  the  better-looking  pieces  were  found  to  exhibit  short 
lifetime  Nd  sites  only.  Some  small  pieces  were  found  which  have  regions 
which  indicate  only  long  lifetime  sites,  but  these  require  excitation  at 
selected  surface  locations.  One  might  conclude  that  either  the  sample 
compositions  are  not  quite  correct  or  that  some  other  means  of  crystal 
growth  must  be  tried. 

In  order  to  gain  further  information,  DTA  data  were  obtained  for 
several  compositions  and  the  thermal  transitions  in  heating  curves  run  at 
10°C/min  are  shown  in  Figure  11.  The  major  transition  in  each  of  the 
compositions  K^LiNdFg,  K^LiNdF^,  K^LiNdFg  and  K2Li2NdF^  occurs  between 
550°  and  575°C.  In  the  samples  of  the  most  interest,  K2LiNdFg  and 
KgLiNdF^,  additional  transitions  are  found  at  760°  and  782°C,  respectively. 
Surprisingly,  the  larger  of  the  thermal  effects  occurs  at  the  lower  tem- 
perature. One  would  have  expected  a melting  point  to  be  more  pronounced 
than,  say,  a phase  transition.  If  the  lower  effect  represents  the  decom- 
position of  a phase,  this  would  suggest  crystallization  of  this  phase 
from  a solvent  as  the  most  reasonable  route  to  crystal  growth.  Certainly, 
the  DTA  data  do  not  suggest  any  of  the  compositions  to  be  congruently 
melting. 
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1.  K2UNdF6 


Figure  11.  DTA  Heating  Curves  for  Several  Compositions 


4. 


DISCUSSION 


The  fact  that  an  unexpectedly  interesting  phase  can  be  made  in  the 
KF-LiF-NdF^  system  is  well  established.  This  is  important  because,  with 
a lifetime  of  280  ps,  it  is  certainly  a member  of  the  small  family  of 
stoichiometric  Nd  materials  in  which  the  active  ion  separation  is  large 
enough  that  quenching  is  much  reduced  compared  to  glasses  and  conventional 
laser  host  crystals. 

A comparison  of  the  emission  spectra  of  highly  crystalline  NdP^Oj^^ 
and  a sample  of  K2LiNdFg  powder  is  shown  in  Figure  12.  The  data  are  taken 
under  essentially  identical  experimental  conditions  and  the  two  materials 
are  seen  to  be  very  similar,  with  the  pentaphosphate  1.05p  emission  some- 
what more  intense.  The  branching  ratio  (simply  defined  as  the  relative 
intensity  of  the  highest  peak  in  the  emission  region  divided  by  that  of 
the  highest  peak  in  the  resonance  region)  of  the  fluoride  is  about  9, 
compared  to  about  4.5  for  the  pentaphosphate.  When  one  considers  that 
the  fluorescence  lifetime  of  the  fluoride  is  more  than  twice  that  of  the 
pentaphosphate  (280  ps  vs.  130  ps),  it  is  found  to  compare  very  well  as 
a potential  mini-laser  material. 
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Figure  12.  A Comparison  of  the  Emission  Spectra  of  NdPgO^^^  and  a Powder 

Seimple  of  the  K^LiNdFc  Composition.  These  represent  ecsentially 
iaentical  experimental  conditions. 
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Further  structural  and  quenching  experiments  were  somewhat  frustrated 

by  the  fact  that  when  dilution  experiments  were  performed,  it  was  found 

that  the  desired  Nd  phase  is  very  restrictive  in  the  rare  earth  ions  it 

will  accept.  For  example,  K2LiNd^Laj^_^Fg  powder  diffraction  data  were 

essentially  those  of  the  tetraf luoride  material,  and  this  conclusion  was 

supported  by  spectroscopic  data.  Gadolinium  and  yttrium  were  tried  as 

diluent  ions  but  the  x-ray  data  showed  obvious  departures  from  the  desired 

phase  with  very  little  substitution  for  neodymium.  Finally,  cerium  was 

found  to  be  much  better  and  samples  diluted  with  Ce^^(K2LiNd  Ce^^ 

appeared  viable  spectroscopically  down  to  about  x = 0.1.  Of  course,  some 

care  must  be  exercised  to  prepare  materials  under  conditions  in  which 
4 + 

Ce  contamination  does  not  occur. 

The  most  interesting  feature  of  Figure  7 in  Section  3,  Results,  which 
shows  the  dependence  of  fluorescence  lifetime  on  concentration,  is  the 
apparent  lack  of  variation  of  fluorescence  lifetime  from  x=0.1tox=l. 
The  material  would  therefore  be  very  interesting  to  study  in  more  detail 
when  single  crystals  and  a crystal  structure  are  obtained.  Of  course,  as 
pointed  out  earlier,  the  compositions  probably  represent  a mixture  of 
KNdF^  and  the  long  lifetime  phase.  If  only  two  Nd-containing  phases  are 
present,  the  data  are  probably  fairly  reliable  since  the  curve  representing 
the  dependence  of  the  fluorescence  lifetime  of  Type  I (the  tetrafluoride) 
is  reasonable  for  a dilution  of  that  material.  It  is  possible,  however, 
that  Nd^^  is  preferentially  found  in  the  long  lifetime  phase.  If  so,  the 
samples  could  be  much  less  dilute  than  indicated.  An  argument  can,  how- 
ever, be  made  that  some  Ce^^  is  incorporated  because  the  branching  ratios 
are  gradually  reduced  from  about  9 for  the  pure  Nd  material  down  to  about 
2 for  X = 0.01.  A regular  dependence  on  Nd  concentration  was  not  found 
for  these  branching  ratios,  but  this  is  not  surprising  because  both  x-ray 
and  spectroscopic  data  indicated  that  the  percentage  of  long  lifetime 
phase  is  probably  not  constant  throughout  the  series.  Still,  the  relative 
increase  upon  dilution  of  the  resonance  emission  in  the  850  nm  to  900  nm 
range  compared  to  that  of  the  emission  in  the  1030  nm  to  1080  nm  range 
is  consistent  with  a reduction  of  self  absorption.  This  indicates  that 
the  sample  is  indeed  diluted  with  Ce  ions.  Thus,  the  conclusion  that 
very  little  quenching  occurs  in  this  material  appears  qualitatively 
correct. 

The  results  of  the  investigation  of  the  interesting  Nd  phase  are 
admittedly  incomplete  at  the  close  of  this  contract.  A continued  investi- 
gation should  be  made  until  some  of  the  questions  presented  here  are 
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cleared  up.  If  a route  to  good  quality  single  crystals  can  be  found,  this 
phase  could  be  especially  valuable  for  elucidating  quenching  mechanisms. 
Its  behavior  is  apparently  quite  different  from  other  stoichiometric  Nd 
compounds’**  and  may  provide  valuable  insight  for  the  development  of 
improved  materials. 
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SECTION  III 

CONCENTRATION  QUENCHING  IN  Nd^'^  STOICHIOMETRIC  MATERIALS* 


♦Reprinted  from  Optics  Communications,  Vol.  23,  No.  3 (December  1977) 
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Basic  experimental  facts  concerning  the  luminescence  of  Nd;(Lai_j(PsO|4  are  reviewed.  Interpretation  of  the  data  by 
a diffusion  model  leads  to  the  possible  presence  of  ~10’’’  cm  ^ energy  sinks  in  this  material.  Consequences  and  weaknesses 
of  this  model  are  discussed. 


The  discovery  that  some  materials  containing  Nd^* 
in  stoichiometric  proportion  (at  concentrations  exceed- 
ing 10-’  ions  per  cm^)  retain  their  luminescence  effi- 
ciency has  given  rise  to  a new  technology  of  miniatur- 
ized, optically  pumped  lasers  [1,2].  At  the  same  time, 
these  developments  have  revived  interest  in  the  mecha- 
nisms of  concentration  quenching  or  the  relative  lack 
of  it  in  materials  such  as  NdP5  0|4. 

It  is  the  purpose  of  the  present  note  to  examine 
which  particular  features  of  the  luminescence  require 
an  explanation,  to  explore  the  possibility  of  providing 
such  an  explanation  in  terms  of  established  models  of 
concentration  quenching,  and  in  general  to  point  out 
the  unsatisfactory  state  of  affairs  concerning  the  under- 
standing of  this  phenomenon. 

Since  most  data  reported  in  the  literature  and  ob- 
tained by  us  pertain  to  Nd^La^_|P50j4,  we  shall  con- 
sider primarily  this  material.  The  following  are  experi- 
mental facts  which  appear  to  dominate  the  lumines- 
cence behavior  of  Nd  ultraphosphate  (NdUP);  other- 
wise also  known  as  pentaphosphate: 

1 ) Concentration  quenching  in  NdUP  is  unusually 
weak.  For  x = 0.1  the  measured  decay  time  does  not 
deviate  from  the  radiative  lifetime  by  more  than  20 
percent  and  at  half  this  concentration  is  essentially 
completely  radiative  [3]. 

2)  The  quenching  constant  varies  linearly  with  the 
molar  fraction  x.  This  fact  has  been  noted  by  Singh  et 
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al.  [4]  and  by  Chinn,  Hong  and  Pierce  [1].  If  we  de- 
note by  r(x)  the  observed  decay  constant  at  the  com- 
position X and  by  fg  the  decay  constant  at  infinite  di- 
lution (x  -►  0),  we  can  define  a quenching  rate  by 

q(.x)=  1/t(x)  - I/tq  ■ 

A normalized  quenching  constant  is  then  defined  as 
4„(x)  = (to-t(x))/t(x). 

Chinn  et  al.  [1]  have  noted  that  for  a whole  family  of 
stoichiometric  materials  (x)  varies  linearly  with  x. 

Our  own  data  obtained  for  just  three  Nd  pentaphosphate 
compositions  confirm  this  and  are  shown  on  fig.  1 . 

The  linear  quenching  appears  to  be  a fundamental  and 
distinguishing  feature  of  these  materials,  yet  some  very 
recent  publications  [5]  do  not  acknowledge  its  exis- 
tence and  attempt  to  explain  quenching  in  terms  of  a 
q^  ~x^  law,  characteristic  ofNdiYAG  [6]. 

3)  The  decay  is  logarithmic  over  most  if  not  all  of 
the  time  accessible  to  measurement.  Decay  curves 
taken  at  the  three  values  of  x and  extending  over  more 
than  two  decades  are  shown  in  fig.  2.  These  we  ob- 
tained by  using  a 5 ns  dye  laser  pulse  exciting  the  sam- 
ples in  the  '*19^2  '*  ^^1/2  transition  and  observing  the 
decay  of  the  ‘*^3^2  "*'^hl/12  luminescence  using  a 
PAR  162  boxcar  integrator  with  a logarithmic  output. 
The  decays  were  followed  to  times  of  the  order  of  1 ms 
after  the  exciting  pulse.  Visual  observation  of  oscillo- 
scope traces  indicated  no  anomalies  at  shorter  times. 
According  to  ref.  [7]  the  decay  at  x = 0.05  is  also  ex- 
ponential. 
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l ig.  1.  Normalized  quenching  constant  versus  molar  fraction 
X in  Nd;(I.ai_;(PsOi4.  (tq  = 320  us;  ref.  [ 1 ]). 

4)  The  decay  time  at  any  concentration  of  Nd  is  in- 
dependent of  temperature. 

Any  model  of  concentration  dependence  of  lumines- 
cence has  to  account  for  all  these  facts  simultaneously. 
We  shall  first  inquire  how  this  can  be  accomplished  by 


using  established  models  of  concentration  quenching. 

Exponential  decay  at  low  concentration  (x  -*  0) 
may  be  due  to  single,  noninteracting  ions.  The  persis- 
tence of  exponential  decay  at  larger  values  of  x implies, 
however,  a mechanism  of  energy  redistribution  which 
maintains  the  equivalence  of  every  excited  ion  in  the 
lattice.  It  is  customary  then  to  treat  such  cases  by  as- 
suming that  the  probability  of  excitation  <i>(r,  r)  is  gov- 
erned by  a diffusion  equation  [8,9] 

d(p(r.  t)/dt  = Dv^0(r,  t) 

- C(/>(r,  r)  S ' — 0(r.  f).  (1) 

n \r-r„\^  ’’0 

In  this  equation  D is  the  diffusion  constant  for  energy 
migration  among  the  Nd  ions,  C is  the  dipole-dipole 
transfer  constant  to  some  unspecified  energy  sinks 
located  at  points  r„ . A solution  of  eq.  (1 ) assumes  a 
simple  exponential  form  if  the  diffusion  is  large.  We 
have  then 

1/t=  I/tq -H/Tp  , (2) 

where 

\ Itq  ^ fi.5  n , (3) 

n is  the  sink  (acceptor)  concentration.  An  expression 
for  the  diffusion  constant  in  terms  of  the  Nd^"^  con- 
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centration  AMs,  following  Forster  [10] 

n = {slTo)R^N*l^ , (4) 

where  s is  a lattice  sum  (of  the  order  of  2.8)  and  Rq  is 
the  Forster  critical  radius  for  Nd - Nd transfer. 
Combining  eqs.  (3)  and  (4)  we  obtain : 

^ = 8.5C'/"(/^4)''\7V.  (5) 

If  the  origin  of  quenching  is  an  interaction  between 
identical  irrns,  eq,  (5)  explains  the  quadratic  depen- 
dence of  q,i  on  N by  simply  assuming  n N.  This  may 
very  well  apply  to  the  case  of  Nd:YAG  in  which  cross 
relaxation  via  the  ‘*1(5/2  level  is  supposed  to  be  fast  [6]. 
At  some  small  level  of  doping  when  eq.  (2)  breaks  down, 
one  should  then  expect  a departure  from  strict  exponen- 
tiality.  We  are  not  aware  of  any  measurements  on 
Nd:YAG  either  supporting  or  contradicting  this  pos- 
sibility. 

If  interaction  between  identical  ions  is  weak,  as  it 
appears  to  be,  the  case  in  Nd  pentaphosphate  (cross 
relaxation  via  ‘*1(5/2  does  not  occur)  [1  l],eq.  (5) can 
formally  account  for  experimental  facts  1 -4  listed 
above.  We  must,  however,  consider « as  a constant, 
independent  of  x,  and  C as  temperature  independent. 

Leaving  for  the  moment  the  reasonableness  of  such 
a model,  we  can  estimate  the  required  concentration 
of  sinks  n by  the  following  approximate  arguments. 

Rq  can  be  estimated  from  its  definition  which  expresses 
the  Nd-Nd  dipole-dipole  transfer  rate  Was  [10] 

W=(l/To)(/?o/«)^  (6) 

where  R is  the  distance  between  neighbors  taken  here 
as  0.6  nm  1 12] . From  the  value  of  W=  2.4  X 10^  s~' 
estimated  by  Liao,  Weber  and  Tofield  [ 1 3] , we  obtain 
Rq  = 2.6  nm.  It  further  follows  from  ref.  [13]  that  the 
diffusion  constant  Z)  * 8 X 10“®  cm^/s  (this,  inciden- 
tally, is  almost  ten  times  smaller  than  the  diffusion  in 
Nd:Y AG  obtained  from  laser  experiments  [14]  and 
casts  some  suspicion  on  either  or  both  of  these  num- 
bers). 

To  estimate  the  Nd  sink  transfer  constant  C,  we 
shall  assume  (see  below)  that  the  same  rate  applies, 
namely  W = C/R^,  which  would  yield  C=  1 .1  X 10~^® 
cm^/s.  Finally,  we  have  Tp  = TQlq„  (x  = 1 ) = 0.006. 

With  these  values  we  find  from  eq.  (5)n  ~ 10*’  sinks/ 
cm^.  Thus,  even  at  low  Nd  content  (x  = 0.01 ),  the  sink 
concentration  amounts  to  no  more  than  one  in  some 
few  hundred  Nd  ions. 
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The  model  would  tend  to  identify  the  sinks  with 
some  crystal  defect  determined  perhaps  by  the  meth- 
od of  preparation  and  constant  over  the  whole  range 
of  Nd/La  proportions.  It  is  conceivable  that  a Nd  ion 
located  at  or  near  such  a defect  could  have  its  energy 
levels  sufficiently  shifted  as  to  allow  for  rapid  pair  re- 
laxation, normally  ruled  out  for  this  material  [11].  If 
this  were  the  case,  the  rate  of  the  final  transfer  step 
could  be  essentially  the  same  as  all  of  the  previous 
transfer  steps  between  “normal"  Nd  ions.  This  would 
justify  our  way  of  estimating  C. 

Another  possibility  would  be  the  presence  of  hy- 
droxyl groups  which  are  known  to  be  effective  quench 
ers  of  Nd^'*’  luminescence  and  whose  concentration 
may  be  determined  by  the  method  of  preparation  and 
be  common  to  all  of  the  crystals  [15].  Although  the 
value  of  C may  be  different,  the  dependence  on  its 
actual  value  is  rather  weak  (see  eq.  (5)).  It  is  a rather 
remarkable  fact  that  as  the  method  of  crystal  growth 
improves,  the  decay  time  at  x = I tends  to  steadily  in- 
crease (chronologically  r = 60  ps,  115  ps,  and  1 35  ps 
- see  refs.  [ 16] , [7]  and  [17],  respectively).  Aside 
from  the  possibility  that  the  defects  are  being  gradual- 
ly eliminated,  it  raises  the  question  whether  the  de- 
pendence of  the  normalized  quenching  constant  on  x 
for  different  stoichiometric  materials  is  universal  as 
suggested  in  ref.  [ 1 ] , or  is  a family  of  close-lying  par- 
allel curves.  The  diffusion  model  would  favor  a family 
of  curves.  Very  careful  measurements  would  be  re- 
quired to  settle  this  point. 

Again,  a consequence  of  the  model  is  that  at  some 
low  value  of  x,  eqs.  (2)  and  (3)  should  break  down  be- 
cause D cannot  be  considered  as  large  any  more.  In 
this  case,  deviations  from  exponential  decay  are  ex- 
pected to  occur,  and  the  initial  portion  of  the  decay 
should  show  the  Forster  dependence  [4, 18].  This 
has  not  been  reported  in  Nd  pentaphosphate  and,  al- 
though precise  data  may  be  lacking,  it  is  not  likely  to 
be  observed  because  at  low  values  of  x the  decay  is 
radiative  anyway.  It  would  appear,  therefore,  that  at 
some  intermediate  concentration  the  radiative  decay 
of  individual  ions  merges  directly  with  a regime  where 
fast  migration  occurs,  the  exponential  decay  being 
maintained  throughout.  An  analytic  treatment  of  the 
qualitative  picture  is  yet  to  be  developed. 

The  main  drawback  of  the  diffusion  model  as  ap- 
plied to  the  pentaphosphate  is  the  arbitrary  nature  of 
the  defects  and  the  postulated  independence  of  n on  x. 
What  are  the  alternatives? 
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To  obtain  a linear  quenching  rate  and  an  exponen- 
tial decay  over  the  whole  range  of  x without  intro- 
ducing defects,  a major  departure  from  the  formulation 
of  the  problem  via  eq.  (1 ) is  indicated.  At  the  same 
time,  a fast  migration  rate  has  to  be  maintained  in  or- 
der to  retain  the  equivalence  of  all  ions.  This  appears 
to  be  the  main  theoretical  problem. 

Without  claiming  rigor  consider  the  following  possi- 
bility. According  to  ref.  [4] , a quenching  interaction 
in  the  pentaphosphate  is  possible  via  the  “'■^li3/2- 
^*9/2  “*'^*15/2  relaxation.  This  should  be  consid- 
erably weaker  than  the  cross  relaxation  in  YAG.  The 
proahility  of  deactivation  by  dipole  dipole  interaction 
can  be  written  in  the  form  of  eq.  (7): 

w = (l/To)(ro/r)6,  (7) 


where  the  average  is  over  the  spatial  distribution  of 
ions,  ff,  is  the  Forster  critical  radius  [ 1 8]  and  r the 
separation  of  Nd  ions.^  long  as  w depends  only  on 
the  average  quantity  1 Ir^,  it  is  immaterial  whether  the 
distribution  is  random  or  uniform;  we  have: 


D 3R^ 


(8) 


and 


/?3 


(9) 


We  thus  obtain  a quenching  rate  linear  in  N,  without 
the  introduction  of  any  defects. 

From  experiments  (fig.  1 ) we  knew  that  vv  = 
q„  (jf  = 1)/Tg  = 6 X 10^  s“* ; this  leads  to  a critical  ra- 
dius of  Tq  = 0.54  nm,  which  is  just  about  equal  to  the 
nearest  Nd-Nd  separation.  The  statistical  fraction  of 
single  ions  as  a function  of  x and  the  number  of  con- 
tacts Z to  other  Nd  sites  is  given  by  [19] 

S = x(l  -xf  . 


For  Z = 8,  rouglily  characteristic  of  NdUP,  and  x = 
0.05,  we  obtain  S = 0.33.  Hence  at  this  concentration, 
66  percent  of  all  Nd  ions  are  singles  and  not  likely  to 
undergo  pair  quenching.  At  jc  = 0.01 , 92  percent  are 
singles  and  hence  decay  radiatively. 

This  simplistic  model  requires  an  energy  migration 
which  is  either  fast  over  the  entire  range  of  x or  varies 
as  a considerably  higher  power  of  N than  indicated  by 


eq.  (5).  In  the  latter  case,  we  can  almost  think  of  a 
fast  diffusion  to  be  rapidly  turned  on  when  x approaches 
and  exceeds  ~0.01 , i.e.,  when  single  ions  decaying  ra- 
diatively begin  to  disappear  from  the  scene  [ 10).  It  is 
possible  that  such  long  range  migration  may  be  supplied 
by  a recently  proposed  hopping  mechanism  [20]  or 
coherent  transfer  [21  ] . 

In  conclusion,  we  believe  that  the  ordinary  diffu- 
sion — defect  model  cannot  he  conclusively  proved  or 
disproved  at  the  present  time.  In  addition,  it  requires 
an  explanation  why  the  Forster  decay  (r’/^-law)  is  not 
observed  at  some  value  of  x.  Measurement  of  the  decay 
for  very  dilute  samples,  and  possible  ways  of  detecting 
or  eliminating  defects  would  shed  light  on  its  validity. 
The  defectless  model,  being  very  incomplete,  requires 
much  further  elaboration.  Direct  measurements  of  the 
diffusion  length  and  high  resolution,  line  narrowing 
experiments  similar  to  those  recently  performed  on 
ruby  [21],  could  shed  light  on  energy  migration  rates 
and  mechanisms. 

Numerous  discussions  with  Drs.  E.J.  Johnson, 

S.J.  Nettel,  B,  McCollum,  L.  Andrews  and  S.R.  Chinn 
are  gratefully  acknowledged. 
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1.  INTRODUCTION 


I 

I 

Pulse  phenomena  in  lasers,  such  as  Q-switching,  transient  spiking  and 
I cavity  dumping  are  usually  treated  using  the  rate  equations.*  In  most 

regions  of  interest,  these  equations  form  a set  of  coupled  monlinear  dif- 
. ferential  equations  and  solutions  are  usually  obtained  using  computer 

I numerical  integration.  Unfortunately,  in  such  calculations,  the  physical 

details  become  obscured;  and  it  is  difficult  to  answer  basic  questions  and 
I compare  different  phenomena  without  large  catalogs  of  numerical  data.  We 

have  examined  in  detail  the  rate  equations  for  a four-level  system  and 
have  found  that  by  sacrificing  rigor  and  piecing  together  solutions  with 
t approximations  appropriate  to  different  levels  of  excitation  we  can  provide 

detailed  insight  into  laser  pulse  phenomena.  The  approximations  are  such 
that  the  resulting  expressions  can  be  trusted  to  give  the  qualitative 
behavior.  In  most  cases  the  resulting  expressions  agree  well  with  numerical 
integration.  However,  for  rigorous  practical  solutions  one  should  resort 
to  computer  numerical  integration  with  the  guidance  provided  by  the 
approximate  analytic  expressions. 

We  compare  various  pulse  modes  such  as  Q-switching,  cavity  dumping  and 
transient  spiking.  We  consider  pulse  narrowing  and  extend  the  equations 
to  include  a system  consisting  of  two  types  of  centers  with  the  same 
energy  level  scheme  but  with  different  lifetimes.  Preliminary  considerations 
appeared  to  indicate  that  such  a two  center  laser  system  might  have  enhanced 
pulse  behavior  and  in  earlier  work  we  came  to  suspect  such  centers  occured  in 
certain  stoichiometric  rare-earth  compounds.  The  two  center  analysis  has 
implications  to  glasses  where  many  different  lifetimes  exist. 
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2.  DETAILS  OF  PHYSICAL  MODEL 


Let  us  consider  a lasing  center  that  can  be  characterized  by  a four- 
level  system. 


Figure  1.  Energy  Level  Diagram  for  a Four-Level  Center 

The  principal  transitions  are  shown  by  arrows  and  appropriate  life- 
times. If  favorable  relations  occur  between  the  magnitudes  of  the  respec- 
tive lifetimes,  the  absorption  constant  involving  the  pumping  transition 
41  is  independent  of  pumping  intensity,  and  the  population  of  level  2 can 
be  neglected.  Under  these  circumstances,  only  t^j^  and  tj2  enter  quanti- 
tatively into  the  dynamics  of  the  laser  output  transition  32,  The  lifetime 
t^j^  is  implicit  in  the  absorption  coefficient,  a;  and  ^32  , if  it  is  limited 
by  radiative  transitions  to  level  2,  corresponds  to  a lifetime  for  spon- 
taneous emission  t.  The  transition  32  is  also  characterized  by  an  Einstein 
coefficient,  per  quantum,  per  mode  for  spontaneous  emission: 


B 


Strv  Avt 


From  this  one  can  define  a quantity  q^: 


tB 


Av 


% 


(1) 


(2) 


where  V is  the  volume  of  the  radiation  field  contained  by  the  resonator. 

The  quantity  q^  is  the  number  of  modes  contained  within  the  line  width,  Av, 

and  will  be  found  to  be  useful  in  describing  the  laser  dynamics.  In  our 

12 

numerical  excimples  we  will  assume  a typical  value  of  q^^  = 1.25  x lo  , 
corresponding  to  a volume  of  1 cm^. 
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The  lasing  centers  are  postulated  to  occur  in  a solid  host  where  they 
can  be  characterized  by  an  absorption  coefficient.  The  laser  medium 
is  in  a resonator  which  is  characterized  by  a cavity  loss  rate  given  by: 


^c  = 


where  r_  corresponds  to  the  coupled  output: 
R 


(3) 


Here  i is  the;  resonator  length  and  R is  the  reflectivity  of  the  output 
mirror.  The  second  term  in  Eg.  (3),  r^,  is  the  equivalent  rate  of  any 
additional  losses. 

(-a  S,  \ 

^ ^ * 's) 

where  is  an  absorption  coefficient  that  accounts  for  absorption  and 
scattering  in  the  cavity  and  f^  is  the  fraction  of  light  scattered  or 
diffracted  at  each  incidence  on  the  two  end  mirrors. 

Note  that  all  losses  are  taken  as  continuous.  This  is  generally  not 
the  case  in  practice.  For  events  of  duration  on  the  order  of  the  transit 
time,  2Z/c  or  less,  the  concept  of  a cavity  loss  rate  can  be  expected  to 
breakdown. 

A useful  unit  of  time  is  given  by  the  cavity  lifetime,  t^,  given  by: 


It  should  be  pointed  out  the.t,  in  general,  the  cavity  lifetime  for 
spontaneous  emission  is  much  shorter  than  that  for  stimulated  emission. 
In  what  follows  we  will  neglect  this  complication. 
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3.  LASER  RATE  EQUATIONS  AND  STEADY  STATE  SOLUTIONS 


The  number  of  excited  centers  in  the  cavity,  n,  under  optical  pumping, 
will  be  governed  by: 

^ = w - VBnq  - Y 

where  we  neglect  nonradiative  transitions  and  w is  the  pumping  rate  which 
is  dependent  upon  the  absorption  coefficient  and  the  incident  light  flux, 
q is  the  number  of  photons  of  a given  mode  present  in  the  cavity.  The 
third  term  is  the  spontaneous  emission  over  the  line  width  and  the  second 
term  is  the  stimulated  emission  per  mode.  Using  Eq.  (2),  Eq.  (6)  becomes: 


The  time  dependence  of  the  number  of  photons  in  the  laser  cavity  is 
given  by: 


^ = BnVq 


+ BnV  - 


<8) 


where  the  first  two  terms  are,  respectively,  the  stimulated  and  the  spon- 
taneous emission  per  mode  of  the  radiation  field,  and  the  last  term  is  the 
rate  at  Jhich  photons  leave  the  resonator*  as  laser  output  and  internal 
'.oases.  Using  Eq . (2).  Eq . (8)  can  be  written: 


(10) 


We  define  two  quantities  whose  significance  will  become  obvious 
below: 


w = r q 
o c ^o 


q T 
^o 


(11) 


With  these  definitions  Eq.  (10)  can  be  written  in  the  convenient  form: 


, r 
dp  _ _c_ 

dt  T 


( 


-3-  + 


^)] 


^ = 
dt 


p (q  + 


3o>  - 


r q 
c ^ 


(12) 


For  cw  laser  operation  one  can  easily  obtain  the  steady  state 
solutions : 


^ ""c  q + q^ 

(13) 

n = WT/(q/q^  + 1) 

1 

(14) 

,\2  ^ 4 ^^1' 

(15) 

In  Fiqure  2 wr  plot  q/q^  versus  w as  given  by  Eq.  (15).  For  w = the 

photon  number  shows  the  virtual  discontinuity  character ist ic  of  the  thresh- 
1 i I f ■ r laser  a " t ion . 
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The  laser  output  intensity  can  be  written; 


if  the  output  coupling  r_  is  optimized: 

K 


For  laser  operation  these  steady-state  solutions  with  optimized  output 
coupling  represent  the  maximum  intensities  available  for  a given  pumping 
rate.  This  limitation  results  from  the  fact  that  in  cw  operation  the  round 
trip  gain  cannot  exceed  the  round  trip  loss.  In  pulsed  operation,  this 
constraint  does  not  exist  and  possibilities  exist  for  obtaining  much  larger 
intensities  than  result  from  Eq.  (18) . Other  applications  require  narrow 
pulses  and  transient  effects  can  give  an  output  pulse  narrower  than  the 
pumping  pulse.  To  understand  these  phenomena  we  examine  approximate  analytic 
solutions  to  the  rate  equations  in  detail  in  the  following  sections. 


4.  BUILDUP  OF  ELECTRONIC  EXCITATION 


The  buildup  of  electronic  excitation  initially  occurs  independent  of 
the  number  of  photons  in  the  cavity.  In  this  situation  Eq.  (7)  can  be 
written ; 


(19) 


This  has  the  general  solution: 


where  C is  a constant  of  integration. 

As  long  as  the  turn-on  time  is  fast  compared  to  t and  w is  reasonably 
flat  during  the  pulse,  little  generality  is  lost  in  taking  w in  the  form 
of  a step  function.  In  this  situation,  the  electronic  excitation  is 
given  by: 


n(t) 


(21) 


where  n{o)  is  the  electronic  excitation  at  t = o. 

There  are  two  aspects  of  the  electronic  excitation  that  are  of  interest 
for  the  laser  dynamics.  There  is  the  total  electronic  excitation  itself, 
i.e.,  n,  which,  when  multiplied  by  the  photon  energy,  represents  the  total 
electronic  energy  available  for  optical  output.  In  addition,  there  is  the 
quantity  p defined  above  which  when  multiplied  by  q gives  the  stimulated 
emission.  Therefore,  p is  the  stimulated  emission  rate  and,  as  such,  plays 
a central  role  in  the  buildup  of  photons  in  the  cavity.  As  we  will  see 
later  p = r^  is  the  threshold  condition  for  exponential  growth  in  the 
photon  population  which  is  required  if  a laser  pulse  is  to  take  place. 
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5.  PHOTON  BUILDUP 

The  presence  of  electronic  excitation  induces  photons  to  build  up  in 
the  cavity.  At  low  levels  the  presence  of  photons  has  virtually  no  .effect 
on  the  electronic  excitation.  At  higher  levels,  coupling  occurs  and  the 
problem  becomes  more  complicated.  In  this  section  we  assume  a variety  of 
time  dependences  for  the  electronic  excitation  and  examine  how  the  photons 
build  up  in  the  cavity. 

For  these  purposes  we  take  the  second  of  Eg.  (12)  : 

^ « (p  - r^)  q + p (22) 

where  p for  present  purposes  is  considered  to  be  a given  function  of  time 
independent  of  the  photon  number,  q.  The  general  solution  to  Eq.  (22)  is 
given  by: 


q 


(/e"^pdt  + cj 


(23) 


where: 

B = /(p  - r^)  dt 

and  C is  a constant  of  integration. 

Solutions  will  differ  depending  upon  the  form  taken  by  p . We  consider 
separately  p a constant,  and  p linear  with  time. 

As  a first  example  consider  p = constant  and  q = q(o)  at  t = 0.  Under 
these  conditions  Eq.  (23)  becomes: 


(P 

q « q (O)  e 


''c> 


1 


(24) 


Thin  •!  '«  If  I'l  lur*'  4 ♦ hf*  » ^ r •*nh' » I 
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Substituting  Eq.  (27)  into  Eq.  (23)  we  obtain: 


The  integral  in  Eq . (23)  can  be  handled  using  the  error  function.  However, 
to  provide  deeper  physical  insight,  we  consider  simple  approximations  suitable 
for  separate  regions.  Schematically  we  show  the  form  the  integral  ta)ces  in 
Figure  5. 


and  the  solution  appropriate  for  t << 


becomes ; 


t / t 
cl  c , 
q = :^  Ve  - 1 
o 


(31) 


which  is  quadratic  for  t <<  t^  and  linear  for  t >>  t^. 


For  t - 


t t 
/ o c 


o si  2 
integral  is: 


)<t<  1 

f 

t 

+ 

l^o^c  ] 

° 

2 / 

a suitable  approximation  to  the 


/ 


e^dt  (t-t)  -|(t-t)^ 

2a  ' o 6 o 


(32) 


and  the  solution  for  q in  this  region  is: 


a 

q ~ /2  t 


1 . *^c  2 ^o  , ^ " ^o  . 1 /tt  ^c  ” ^o^ 

_l  + _e  + — ^ ^2  2t ;2 

c o c V o t 


(33) 


/t  t 

Finally  for  [ t^  + I <<  t 


/o-“ 


(34) 


giving 


q ♦ 1 


-1  !2\  1 

,,  ^o  , ^ ^c\  2" 

t * t " * ** 

l V r o 


^o>' 


(35) 
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where  p can  be  a function  of  time.  In  addition,  for  t >>  t^  one  could 
expect  the  right  side  of  Eq.  (36)  to  approach  a good  approximation  to  q. 
Indeed,  one  can  obtain  such  a solution  by  setting  the  derivative  equal  to 
zero  in  Eq.  (22) . We  assume  that  in  some  range  of  t,  Eq.  (35)  is  a valid 
solution.  If  p = at,  Eq.  (36)  becomes: 


i-f 

o 


This  agrees  well  with  Eq.  (31)  for  t^  <<  t <<  t^,  and  can  be  expected 
to  be  a good  approximation  at  higher  values  of  t until  dq/dt  increases  too 
rapidly  as  t approaches  t^. 

Now  we  consider  how  these  solutions  might  be  applied  to  the  results  of 
the  last  section.  There  we  found  that  if  we  apply  a step  function  pumping 
pulse,  the  electronic  excitation  can  be  expected  to  follow  Eq.  (21)  which 
can  be  written: 


(1  - e-^/^) 


The  straightforward  approach  to  obtain  the  time  dependence  of  q would 
be  to  substitute  Eq.  (38)  into  Eq.  (23).  Unfortunately,  in  doing  so  we 
obtain  integrals  which  are  not  easily  solved.  One  way  to  overcome  this 
difficutly  is  to  write  Eq.  (38)  in  terms  of  a series  giving: 


w , , , n+1  1 /t\* 

■>  ■ 5-2^  '-11  7m  It) 


or  1 f t • ■ T 


If  <<  T only  the  first  term  of  the  series  in  Eq.  (39)  need  be 

considered,  and  all  of  the  expressions  based  upon  a linear  behavior  of  p 

w 

can  be  taken  over  by  setting  t^  = -^  t . 

The  case  where  w < is  also  a simple  situation.  Here  t^  does  not 

exist  since  p does  not  reach  the  level  of  r . Eq.  (31)  holds  for  t < t, 

w c 

nevertheless,  if  t^  is  replaced  by  t giving: 


At  longer  times  this  connects  smoothly  with  an  expression  obtained  from 
Eq.  (36)  with  substitution  of  Eq . (38): 


q 


^ (1 
w 

o 


1 - (1  - e"^/^) 

w 

o 


(42) 


In  the  case  where  w w^  the  situation  is  more  complicated;  but,  as 

it  turns  out,  not  much.  Eqs.  (41)  and  (42)  are  still  valid,  except  that 

the  behavior  given  by  Eq . (42)  is  terminated  abruptly  when  the  stimulated 

emission  rate  reaches  the  cavity  loss  rate  at  t = t . To  obtain  t one 

00 

needs  the  more  general  expression: 


(43) 


Things  happen  so  rapidly  in  the  neighborhood  of  t that  we  can  still 
use  a linear  approximation  for  the  stimulated  emission  rate  where  a is 
given  by: 
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I 

( 


This  expression  is  qualitatively  no  different  than  Eq.  (35).  However, 
the  constants  a and  t^  have  the  special  definitions  given  above  and  the 
constant  b cannot  in  general  be  obtained  in  a simple  manner  from  the 
behavior  at  earlier  times.  If  w , B can  be  estimated  from  Eq . (35): 


(46) 
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NUMERICAL  INTEGRATION  RESULTS  FOR  PHOTON  BUILDUP 


Computer  numerical  integration  was  used  to  obtain  the  time  dependence 
of  g and  these  results  were  compared  with  the  approximate  analytic  expres- 
sions obtained  in  Section  5.  The  numerical  integration  employed  similar 
techniques  to  those  reported  in  several  papers.' 

In  Figure  6 numerical  results  for  w/w  = 0.1,  1.0  and  10  are  shown. 

-21  ° 

At  levels  well  below  g 10  q^,  the  characteristics  of  Eq.  (31)  are 

displayed.  Below  t = t^,  q varies  quadratically  passing  over  to  a linear 
dependence.  In  fact,  the  numerical  results  are  in  exact  agreement  with  the 
analytic  expression. 

At  high  levels  of  q there  is  a region  where  the  cavity  photon  number 
saturates  for  w < w . For  w = w , q takes  on  a superlinear  dependence  in 


this  region  in  accordance  with  that  expected  when  p ^ r^ . An  even  greater 
superlinear  region  occurs  for  w > w^,  which  is  apparently  related  to  the 
onset  of  exponential  behavior  expected  when  the  stimulated  emission  rate 
exceeds  the  cavity  loss  rate. 

-14 

The  region  above  q = 10  is  shown  in  a larger  scale  in  Figure  7. 
Comparison  is  made  between  the  numerical  integration  and  points  calculated 
according  to  Eq.  (31),  which  is  exact  at  lower  levels,  and  calculations 
based  upon  the  asymptotic  Eq.  (36). 

For  w/w  = 0.1  departures  from  Eq.  (31)  are  observed  even  at 
-14  ° 

q = 10  q^.  However,  good  agreement  is  obtained  between  Eq.  (36)and 

the  numerical  integration.  This  agreement  with  the  asymptotic  expression 
should  continue  indefinitely  since  p cannot  exceed  0.1  r^ . 

A similar  situation  occurs  for  w/w^  = 1.0.  However,  since  p approaches 
r^  the  denominator  in  Eq.  (36)  ma)ces  itself  felt  and  ^ becomes  too  large 
to  permit  q to  keep  up  with  Eq.  (36) . However,  since  p cannot  exceed  the 
cavity  loss  rate,  q does  saturate  at  a value  equal  to  /q^  (not  shown  in 
Figure  7 ) . 

For  w/w^  =10,  in  addition  to  the  above  features,  q shows  a region 
where  its  variation  appears  nearly  vertical.  This  is  the  region  where  the 
stimulated  emission  rate  exceeds  r^,  giving  exponential  growth.  This 
region  is  of  principle  interest  for  lasers  and  is  examined  in  detail  below. 

In  Figure  8 we  show  numerical  results  near  t ~ t^,  i.e.,  for  the 
region  where  p % r^.  The  values  of  q at  t = t^  are  found  to  agree  quite 
well  with  those  predicted  by  Eq.  (33),  as  shown  in  Table  1. 


TABLE  1 

COMPARISON  OF  TYPICAL  VALUES  OF 
PHOTON  NUMBER  AT  THRESHOLD 


W/Wo 

q(t^)  FROM  EQ.  (33) 

q(to)  FROM 

NUMERICAL  INTEGRATION 

10 

39.6 

44 

5 

58.2 

64 

2 

103.3 

123 

In  Figure  9 the  various  curves  of  Figure  8 have  been  shifted  vertically 
to  facilitate  comparison  of  their  slopes  at  t = t^.  As  shown  in  the  figure 
the  slopes  are  all  close  to  a value  equal  to  r^,  as  predicted  by  Eq.  (33). 


I 
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Figure  9.  Curves  of  Figure  8 Shifted  to  Compare  Slopes  at  t = t 
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In  Figure  10,  the  nvunerical  integration  shows  clearly  the  transition 

to  exponential  growth  in  the  range  t ss  t . Above  the  transition  q grows 
2 ° 

as  exp  (t  - t^)  in  agreement  with  Eg.  (35)  and  has  an  intercept  at  t = t^ 
of  70,  compared  to  79  calculated  according  to  Eg.  (35). 


Figure  10.  Transition  to  Exponential  Growth  in  Photon  Number 

In  Figure  11  we  show  the  exponential  growth  for  a variety  of  examples 

in  the  range  from  w/w^  = 1.1  to  w/w^  = 10.  It  is  seen  that  Eg.  (45) 

is  a good  approximation  to  the  behavior  with  the  constant  B varying  from 

350  for  w = 1.1  w to  44  for  w = 10  w . 

o o 

The  numerical  integration  shows  that  the  approximate  analytic  expres- 
sions developed  in  Section  5 are  valid  analytic  approximations  for  the 
growth  of  q during  a step  function  pumping  pulse  and  therefore  should 
provide  a good  foundation  for  analyzing  pulsed  lasers. 
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7.  COUPLING  BETWEEN  THE  PHOTON  FIELD  AND 
THE  ELECTRONIC  EXCITATION 


We  have  seen  that  the  first  requirement  for  the  initiation  of  a pulse 
is  that  the  spontaneous  emission  rate  rises  to  equal  the  cavity  loss  rate. 
Under  these  conditions  one  gets  exponential  growth  in  q.  But  how  high  does 
q go  before  it  reaches  a peak  value?  From  the  equations  it  is  apparent 
that  q will  rise  until  the  stimulated  emission  rate  falls  back  to  again 
equal  the  cavity  loss  rate.  Therefore,  the  second  condition  for  a laser 
pulse  is  that  the  stimulated  emission  rate  spend  some  time  where  it  is 
greater  than  the  cavity  loss  rate.  Further,  one  might  expect  that  the 
peak  intensity  of  the  laser  pulse  might  depend  upon  the  peak  level  attained 
by  the  stimulated  emission  rate. 

Let  us  assume  that  by  some  process  the  stimulated  emission  rate  has 
risen  and  leveled  off  at  some  maximum  value,  p . Further  assume  that  q, 

ITldX 

though  greater  than  q^,  is  at  a level  small  compared  to  what  it  will 
attain  in  the  pulse. 

In  this  situation  Eq.  (12)  can  be  approximated  by: 


Ip  = 

dt 


^o 


^ = 
dt 


q 


(p 


(47) 


where,  since  exponential  growth  in  q is  occurring,  we  neglect  the  pumping 
term.  We  can  eliminate  time  from  these  equations  giving: 


^ = _ 
dp 


(48) 


The  solution  to  this  equation  is  given  by: 


= ^^o  ( 


, max 

p “ p - r In  

^max  ^ c p 


(49) 


We  can  immediately  write  for  the  maximum  value  attained  by  qt 


q = q X 

^max 


p - r I 1 + In 
"^max  c 


(i . 


(50) 


where  we  neglect  the  value  of  q at  the  beginning  of  the  exponential  growth. 
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To  obtain  the  time  dependence  of  q we  could  substitute  Eq.  (49)  into 
Eq.  (47).  Since  this  gives  difficult  integration,  we  ma)ce  use  of  the 
approximation : 


q 


max 


1 q T 

1 ^o 

2 r 


c 


(P 


(51) 


This  approximation  is  valid  for  p 
value  of  q. 

Substituting  Eq.  (51)  in  Eq.  (47) 


% r and  therefore  near  the  pea)< 
c 

gives : 


d (p 


q 

^max 

% 


(52) 


Near  the  peak  the  variation  in  (p  - r^)  affects  the  integral  much  greater 
than  the  variation  in  p.  We  set  the  ratio  % 1.  Integration  and 

substitution  in  Eq.  (49)  then  yields: 

q q sech^  ( — | (53) 

^ ^max  \y } 

where  y corresponds  to  the  pulse  half -width  and  is  given  by: 


1 

Y 


1 

2 


Substituting  from  Eg.  (50): 


(54) 


1 

Y 


p 

max 


r 

c 


|l  + In 


^max  j 


(55) 


If  p is  too  large,  one  can  expect  errors  in  Eq.  (53)  near  the  half  power 
indx 

points.  These  errors  will  tend  to  compensate,  making  the  rise  time  sharper 
and  the  decay  time  slower. 


In  the  case  where  p < 2r  , Eq.  (55)  becomes  simplified  to: 
mdx  c 


y 


2t 


c 


r 


c 


p 

^max 


r 

c 


(56) 
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8.  CAVITY  DUMPING 


The  first  pulsed  laser  phenomena  we  consider  is  cavity  dumping.  This 
technique  can  produce  peak  powers  enhanced  over  the  optimum  cw  powers  and 
provides  a means  of  modulating  the  laser  output  in  the  megahertz  range. 

In  cavity  dumping  the  output  mirror  is  replaced  by  a 100%  reflecting 
mirror.  The  circulating  power  is  allowed  to  build  up  and  means  are  pro- 
vided to  quickly  deviate  all  or  part  of  the  circulating  beam  out  of  the 
cavity.  In  terms  of  the  rate  equations,  this  amounts  to  modulating  the 
output  coupling. 

The  number  of  photons  in  the  cavity  corresponds  to  the  cw  stored  energy 
and  is  given  by  Eq.  (16).  Well  above  threshold  and  with  the  output  coupling 
reduced  to  zero,  q is  limited  by  scattering  losses  and  is  given  by: 


q = (wtj  - q^) 


(57) 


The  electronic  excitation  is  given  by: 


n 


T 


(58) 


The  corresponding  stimulated  emission  rate  is: 


(59) 


We  have  a situation  where  there  can  be  much  energy  stored  in  the  photon 
field  in  the  cavity  (if  t^^  is  long) , and  little  stored  in  the  electronic 
excitation  if  (t^^  ^ i)*  if  riow  the  output  coupling  is  abruptly  switched  to 
give  a short  cavity  lifetime,  t^ , a peak  output  intensity  results,  given  by: 


(60) 


The  output  can  result  in  enhanced  intensity  in  two  respects.  First, 
the  circulating  power  is  higher  than  cw;  and  second,  the  output  coupling 
can  be  much  smaller  than  the  optimum  coupling  in  the  cw  case. 


101 


The  decay  is  given  by: 


^ = 
dt 


P 


(61) 


Since  >>  p , the  decay  of  the  output  is  exponential  and  given  by: 


wt . - q t 

I C 


(62) 


Once  the  output  coupling  is  switched  back  to  zero,  the  circulating 


power  recovers  in  a time  t^,  given  by  Eg.  (43) : 


wt , 


t = tin  — 

o wt.  - 


(63) 


Had  the  laser  been  operating  cw  at  some  optimum  output  coupling,  t^,  it 
would  have  had  an  output  given  by  Eq.  (18).  The  optimized  cavity  dumped 
output  is  given  by: 

^ (64) 

Therefore,  cavity  dumping  increases  the  peak  power  by  the  factor: 


F = 


(65) 


In  principle,  an  output  coupling  of  100%  would  be  possible  giving  t^ 
equal  to  the  round-trip  time  of  the  cavity.  However,  practical  considera- 
tions in  real  systems  limit  the  maximum  output  coupling  to  50%  as  used  in 
Eq.  (62) . With  high  output  coupling  the  decay  departs  from  simple 
exponential  behavior  showing  the  effect  of  discrete  passes  in  the  resonator. 
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9.  Q-SWITCHING 


Q-switching  is  another  pulsed  laser  phenomenon  where  one  modulates 
the  cavity  loss  rate.  In  this  case  one  initially  makes  the  cavity  loss 
rate  very  large  by  introducing  losses.  This  permits  the  stimulated  emis- 
sion rate  to  rise  to  high  levels  without  lasing.  However,  since  the  losses 
are  high  the  number  of  photons  in  the  cavity  are  low.  This  situation  is 
exactly  the  opposite  of  that  in  the  cavity-dumping  situation.  If  we  now 
switch  the  cavity  loss  rate  to  a value  much  less  than  the  stimulated 
emission  rate,  exponential  growth  in  q occurs  resulting  in  an  intense 
output  pulse. 


To  analyze  this  process  we  suppose  a step  function  pumping  pulse  is 

applied  to  the  laser.  Since  initially  the  losses  are  high,  the  laser  does 

not  reach  threshold.  The  stimulated  emission  rate  rises  according  to 

Eq.  (38)  , therefore  reaching  a maximum  of  w/q^  in  times  comparable  to  a 

radiative  lifetime  t.  Since  p never  exceeds  the  cavity  loss  rate,  the 

photon  number  saturates  at  the  value  — t . To  the  extent  to  which  t can 

^o  c 

be  considered  to  approach  zero,  q can  be  considered  zero.  The  electronic 
excitation  saturates  at  a value  of  wt. 

If  now  r^  is  switched  to  a value  much  less  than  w/q^,  the  photon 
number  rises  exponentially  according  to  Eq.  (24) : 


q = 


r (p  - r )t 

-£ e - 1 

- r 

c L J 


(66) 


This  exponential  rise  in  q will  continue  until  q becomes  comparable  to  the 
saturated  electronic  excitation,  and  p starts  to  become  depleted.  Although 
q may  be  less  than  q^  it  quickly  exceeds  q^,  so  that  we  can  take  over  the 
results  of  Section  7.  From  Eq.  (50)  we  can  calculate  the  maximum  value 
attained  by  q; 


(67) 


q 5=  WT 

w 

1 - ^ 

/l  + In 
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The  peak  output  is; 
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w \ 

" ^)\ 

(68) 


A lower  limit  to  the  pulse  width  is  obtained  from  Eq.  (54) : 


As  pointed  out  previously,  the  pulse  half-width  given  by  Eq.  (54)  is  based 

on  a symmetric  pulse  which  is  a good  approximation  only  for  low  values  of 

excitation.  The  situation  for  Q-switching  is  li)cely  to  be  more  like  that 

shown  schematically  in  Figure  13,  with  a steep  rise  and  a slow  decay 

-r  t 

eventually  going  as  e ^ . This  consideration  would  suggest  that  the 
pulse  width  would  be  larger  than  given  by  Eq.  (54)  with  an  upper  limit: 

Y ~ (72) 

This  value  of  half-width  is  consistent  with  a pulse  energy  given  by: 

E = WT  (73) 


t/t 


c 


Figure  14.  Numerical  Integration  Results  for  Q-Switch  Pulse 
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In  Figure  14  we  show  the  results  of  numerical  integration  for  the 
q-switch  situation.  The  general  conclusions  obtained  above  are  confirmed. 
The  pulse  is  asymmetric,  as  expected,  with  a half-width  t^.  The  photon 
number,  q,  rises  to  a value  comparable  to  the  electronic  excitation.  A 
calculation  of  q„^^  using  Eg.  (67)  gives  q„,„  = 5.92  x lo^  q which  agrees 

niaX  in«X  O 

well  with  the  numerical  integration. 


Figure  15.  Variation  of  q with  Spontaneous  Lifetime 
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10.  TRANSIENT  SPIKE  LASER  OPERATION 


We  now  consider  the  transient  spike  that  occurs  at  the  initiation  ol 
optical  pumping.  In  this  form  of  pulsed  operation  we  modulate  the  pumping 
rate,  w.  Our  discussion  will  be  limited  to  the  application  of  a step  func- 
tion pumping  pulse.  Such  a situation  applies  to  the  practical  situation 
where  the  pumping  pulse  is  reasonably  flat  and  where  the  turn  on  time  is 
short  relative  to  the  spontaneous  lifetime. 

As  in  other  forms  of  laser  pulses,  a pulse  is  initiated  when  the 
stimulated  emission  rate  rises  to  the  order  of  the  cavity  loss  rate.  There 
are  three  events  during  the  optical  pulse  that  have  particular  significance 
First,  there  is  the  point  in  time  where  p first  exceeds  r^.  This  situation 
permits  exponential  growth  in  q.  However,  since  q starts  from  a very  low 
level,  the  stimulated  emission  rate  is  permitted  to  rise  significantly 
above  r^.  There  is  little  laser  output  during  this  period.  Second,  there 
is  the  point  where  p reaches  its  maximum.  The  photon  number  is  still  at  a 
relatively  low  level,  but  the  maximum  value  achieved  by  the  stimulated 
emission  rate  determines  the  ultimate  strength  of  the  pulse.  Third,  there 
is  the  point  at  which  p begins  to  drop  below  r^.  Coincident  with  this,  q 
reaches  its  peak  and  the  output  pulse  begins  to  decay  in  intensity. 

The  stimulated  emission  rate  during  the  initial  period  is  given  by 
Eq.  (38).  However,  in  general,  the  pulse  is  so  short  that  one  can  use  the 
linear  approximation  with  the  constant  a given  by  Eq.  (44)  . 


(74) 


and  q during  the  initial  period  is  given  by  Eq.  (45) : 


q + 1 = B exp 


(75) 


where  B depends  upon  w/w^  and  can  be  determined  as  discussed  in  .Section  5. 

The  condition  that  p reaches  a maximum  is  that  ^ = 0 in  Eq.  (12) 
which  is  satisfied  by: 


(76) 
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The  maximum  stimulated  emission  rate  is  given  by  the  simultaneous 
solution  of  Eqs.  (74),  (75)  and  (76).  This  solution  is  given  approximately  by: 


(p  - r ) ■ 
max  c 


(»^  - 


2 - 


In 


q ( J!L  _ i\ 

^O  I p w / 

\ max  o / 


(77) 


If  w >>  V7^,  in  accordance  with  Eq.  (46)  the  log  function  can  be 
replaced  in  this  and  subsequent  expressions  by: 
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(78) 


To  obtain  high  stimulated  emission  rates  one  should  use  a high  pumping 
rate  and  select  material  with  a low  spontaneous  lifetime.  However,  as  t 
tends  to  zero  the  stimulated  emission  rate  has  a limit  given  by: 


limit 


T -►  0 


p w w 

^max  = — r = — 
w c q 
o ^o 


(79) 


This  limit  is  approached  closely  for 


T < 


In 


("'o  ' 


^ - 1 
w o 

o ^ 2'n 


(80) 


w 

w 


by: 


For  higher  values  of  r where  p r , Eq,  (77)  can  be  approximated 

max  c 


P = r + r \ 2 — 
max  c c 


? it  - 
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In  the  range  of  low  spontaneous  lifetimes  specified  by  Eq.  (80),  the 
peak  value  of  q is  obtained  by  substitution  in  Eq.  (50)  is  given  by: 


q 

^max 

T - 0 


1 + In 


(82) 


In  this  range  we  see  that  increases  linearly  with  lifetime.  At 

longer  lifetimes  approaches  r^ , and  Eq.  (50)  can  be  approximated  by: 


q It 
o c 


"■max 


(P 
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Substitution  of  Eq.  (81)  gives: 


(83) 


(84) 


(85) 


If  w >>  w we  have  from  Eq.  (46) : 
o 


(86) 


Expressions  (85)  and  (86)  provide  the  principal  dependence  of  q„^„  on 
lifetime  in  the  range  where  Eq.  (84)  is  valid.  In  this  range  the  dependence 
is  very  weak  relative  to  the  linear  dependence  at  lower  values  of  t.  At 
intermediate  values  of  t Eq.  (85)  dominates  the  behavior,  continuing  the 
trend  toward  greater  values  of  q . Eventually  this  trend  is  reversed 
by  the  variation  of  B with  lifetime,  which  dominates  the  behavior  at  high 
values  of  lifetime.  This  behavior  is  shown  graphically  in  Figure  15. 


Ill 


t/1.  AFTER  TRIGGER 


Figure  16.  Numerical  Integration  for  Optical  Spiking 

We  now  consider  other  aspects  of  the  pulses  in  the  range  of  the  maximum 

in  Figure  15  where  we  can  assume  p r . The  maximum  laser  output 

max  c 

intensity  becomes; 


With  optimized  coupling  this  becomes: 


The  enhancement  factor  of  the  peak  intensity  over  cw  is: 


From  Eq.  (53)  the  time  dependence  of  the  pulse  is  given  by: 


(87) 


(88) 


(89) 
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where  the  pulse  half-width  is  given  by: 


y = 2.5 


and  a pulse  energy  given  by: 


E = 2 . 5 /q  r T 
J^o  c 
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(92) 


We  can  check  these  conclusions  in  another  way  by  integrating  Eq.  (47) 
term  by  term.  From  the  first  of  Eq.  (47)  we  obtain: 

iq^  f ^ f - y* qpdt  - J pdt  (93) 

We  take  as  the  initial  limit  of  integration  the  point  at  which  P = p. 


and  the  second  limit  at  p 


max 


The  term  on  the  left  becomes: 


iq  / dp  = Tq  (p  - r ) (94) 

^oj  ^o  ^max  c 

For  a short  pulse  the  first  term  on  the  right  of  Eq.  (93)  is  negligible. 

We  assume  that  during  the  pulse  p does  not  deviate  much  from  r^.  The 
second  term  becomes : 

/r^qdt  = | (95) 


where  E is  the  energy  in  the  pulse  and  we  assume  for  simplicity  that  the 
pulse  is  symmetrical.  We  further  neglect  the  last  term,  assuming  q becomes 
somewhat  larger  than  q^  during  the  pulse.  Therefore,  we  have: 


E = 2Tq  (p  - r ) 
max  c 


(96) 


113 


Substituting  from  Eq.  (81): 


This  is  in  substantial  agreement  with  Eq.  (92). 

The  analysis  in  this  section  made  liberal  use  of  the  results  obtained 
in  Section  7 which  were  based  upon  the  assumption  that  at  the  time  the 
stimulated  emission  rate  reached  its  maximum,  q was  somewhat  greater  than 
q^  but  at  a smaller  level  compared  to  what  it  would  attain  in  the  pulse. 

To  test  the  self-consistency  of  these  assumptions,  we  calculate  the  values 
of  these  quantities  for  several  examples  in  Table  2. 

TABLE  2 

TYPICAL  VALUES  OF  VARIOUS  QUANTITIES  DURING  TRANSIENT 
SPIKE  OPERATION  FOR  T = 10^  t 

c 


w/w 

o 

r p 
c max 

q/q  at  p = p 

n/  'Iq  ^max 

F 

10 

70 

1.22 

7.2 

25.8 

5 

120 

1.14 

3.38 

24.5 

2 

190 

1.07 

0.87 

22.6 

1.5 

200 

1.05 

0.42 

21.8 

1.1 

450 

1.04 

0.06 

19.4 

For  the  assumptions  to  be  satisfied,  the  numbers  that  appear  in  Column  4 

should  be  greater  than  unity  and  small  relative  to  those  in  Column  5.  The 

first  condition  is  satisfied  for  w > 2w  and  the  second  in  all  cases.  There- 

o 

fore,  the  analytic  expressions  obtained  should  be  good  approximations  except 

for  w too  close  to  w . 

o 

In  Figure  16  we  show  numerical  integration  results  for  transient  spi)<es 

for  two  values  of  spontaneous  lifetime  differing  by  an  order  of  magnitude. 

2 

The  pulses  are  reasonably  symmetrical  as  expected  for  a sech  yt  and  the 
pealts  are  nearly  independent  of  lifetime,  as  shown  in  Table  3.  Values  of 
*^max  half-width  agree  quite  well  with  the  corresponding  values  calculated 
from  the  approximate  analytic  expressions.  In  Table  3 comparison  is  made  of 
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values  of  and  y calculated  from  analytic  expressions  and  values  obtained 

'*ITiaX 

through  numerical  integration.  Reasonable  agreement  occurs. 


TABLE  3 

COMPARISON  OF  ANALYTIC  EXPRESSIONS  AND  NUMERICAL 
INTEGRATION  FOR  OPTICAL  SPIKE  OPERATION 


3 

T = 10  t 

c 

T = 10^  t 

c 

q /q  • Calculated 

^max  ^ , 

• Numerical 

Integration 

182 

192 

213 

238 

y/t  • Calculated 

16.4 

51.9 

• Numerical 

Integration 

16 

52 
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11.  PULSE  NARROWING 


Frequently  the  situation  arises  where  one  has  a laser  output  that  has 
more  than  enough  intensity  for  the  application.  However,  it  may  not  be  the 
right  wavelength  and  the  pulse  width  may  be  excessive  for  the  application. 
Such  may  be  the  case  with  a Q-switched  glass  or  ruby  laser,  or  with  a molec- 
ular nitrogen  laser.  One  can  use  such  a laser  to  pump  a second  laser  whose 
output  occurs  at  a longer  wavelength  and  where  interesting  pulse  narrowing 
phenomena  may  occur. 

Of  course,  any  of  the  pulse  phenomena  we  have  talked  about  is  capable 
of  pulse  narrowing,  since  for  any  laser  pulse  that  exists  there  is  a pumping 
pulse  that  is  longer.  As  we  have  seen,  cavity  dumping  and  Q-switching  gives 
pulse  widths  on  the  order  of  a cavity  lifetime.  However,  the  spontaneous 
lifetime  may  be  too  short  for  Q-switching  or  for  some  other  reason  one  can 
not  or  would  prefer  not  to  use  such  techniques.  In  this  case,  one  can 
exploit  the  properties  of  transient  spike  operation. 

The  use  of  the  transient  spike  approach  to  pulse  narrowing  was  first 
suggested  by  Roess^  who  demonstrated  the  effect  by  numerical  integration 
of  the  rate  equations.  The  technique  was  further  developed,  largely  in  an 
emperical  manner,  by  others. 


In  our  treatment  of  transient  spike  operation  for  long  spontaneous 
lifetimes,  we  obtained  the  pulse  width  given  by  Eg.  (91). 


(98) 


Therefore,  for  narrow  pulses  one  should  use  short  cavities  with  high 
output  coupling  and  pump  as  high  above  threshold  as  possible. 


If  this  pulse  width  is  to  be  isolated  from  additional  oscillations  or 
significant  cw  output,  the  pumping  pulse  must  be  turned  off  in  a time  slightly 
longer  than  the  delay  that  occurs  between  the  initiation  of  pumping  and  the 
transient  spike.  This  delay  is  given  by  Eq.  (43) : 


t 


o 


(99) 


Returning  to  Eq.  (98)  we  note  that  for  narrow  pulses  it  is  helpful  if 
the  spontaneous  lifetime  is  short  as  in  dye  lasers.  However,  if  the  life- 
time approaches  the  cavity  lifetime  the  assumptions  that  led  to  Eq.  (98) 


\ 
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are  no  longer  valid,  and  we  need  to  derive  a new  expression  to  replace  it. 
For  short  cavity  lifetimes  the  maximum  stimulated  emission  rate 
attained  in  a transient  guise  is  given  by  Eg.  (79) . This  leads  to  a 
maximum  photon  number  given  by: 


g = T 
^max 


w In  — 
o w 


(100) 


From  Eg.  (54)  we  obtain  the  pulse  width  for  short  spontaneous  lifetimes 


Y 


/2  t 


c 


1 


(101) 


For  pumping  3X  threshold,  this  gives  a pulse  width  approximately  egual 
to  the  cavity  lifetime.  This  is  precisely  what  has  been  observed  in  dye 
lasers  where  pulse  widths  less  than  a nanosecond  have  been  observed. 
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12.  TWO  CENTER  LASER 


It  is  interesting  to  speculate  whether  a pulsed  laser  might  profit 
from  two  types  of  lasing  centers  with  the  same  energy  level  structure  but 
with  greatly  different  lifetimes.  One  would  attempt  to  store  large  amounts 
of  energy  in  long  lifetime  centers  and  to  utilize  the  properties  of  the 
fast  centers  to  dump  this  energy  quickly  in  a sharp  pulse.  Perhaps  in  this 
way  one  could  obtain  transient  spike  operation  to  rival  Q-switch  operation 
without  the  complications  of  extra  intracavity  elements.  The  feasibility 
of  this  approach  is  not  obvious  and  it  is  interesting  to  go  through  the 
arguments,  which  are  at  first  encouraging,  but  eventually  show  that  there 
is  no  advantage  to  such  a two-center  approach. 

First  we  consider  how  to  modify  the  rate  equations  to  take  into  account 
two  sets  of  centers  of  spontaneous  lifetimes  t^  and  tj^,  respectively.  We 
assume  tj^  is  something  like  a factor  of  five  or  ten  greater  than  t^. 

If  the  line  width  is  independent  of  the  spontaneous  lifetime  (as  in 
the  case  of  the  Doppler  limited  situation)  and  is  the  same  for  both  types 
of  centers,  then: 
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dt 

= P (q 

1 - 

where  subscripts  a and  b designate  quantities  associated  with  the  fast  and 
slow  centers  respectively,  and  P“Pa'''Pj3- 
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In  a single  center,  laser  p and  n are  tied  together  by  a single 

spontaneous  lifetime.  In  a two-center  laser  where  t and  t,  are  vastly 

a D 

different,  an  additional  degree  of  freedom  exists  where  the  energy  storage 
can  be  dominated  by  the  slow  centers  nj^,  and  p is  dominated  by  the  fast 
centers  through  n /t  even  though  n may  be  small.  It  is  this  separation 
of  roles  in  pulsed  laser  operation  that  we  would  like  to  exploit.  If  one 
is  to  hope  to  take  advantage  of  this  duality,  one  must  consider  the  relative 
magnitudes  of  the  pumping  rates.  The  effect  of  the  relative  magnitudes  of 
the  pumping  rates  on  the  electronic  excitation  is  illustrated  in  Figure  17 
where  the  contributions  of  each  center  to  the  electronic  energy  storage 
are  given  by  solid  lines  and  the  contributions  to  gain  are  given  by  dashed 
lines.  The  first  set  of  curves  are  given  for  w^  > 5Wj^  where  we  assume  a 
ratio  of  5 between  the  lifetimes.  In  this  case  the  laser  dynamics  are 
completely  dominated  by  the  fast  centers  since  the  contributions  of  the 
slow  centers  to  both  energy  storage  and  to  stimulated  emission  rate  are 
small  relative  to  the  contributions  of  the  fast  centers. 

For  w,  < w < 5w,  , we  arrive  at  the  desirable  situation  where  the 
b a b 

energy  storage  is  dominated  by  the  slow  centers  and  the  stimulated  emission 
is  dominated  by  the  fast  centers.  It  is  important  not  only  that  the  stimu- 
lated emission  rate  attain  high  values  but  also  that  it  react  rapidly  to  a 
trigger  pulse. 

We  have  shown  that  one  can  arrange  a situation  where  the  energy  storage 
in  an  initial  below  threshold  pulse  is  dominated  by  the  slow  centers.  In 
a subsequent  intense  narrow  pulse  the  stimulated  emission  rate  can  be  domi- 
nated by  the  fast  centers.  The  question  that  remains,  however,  is  whether 
the  energy  stored  in  the  slow  centers  would  contribute  significantly  to  the 
laser  spike  triggered  by  the  later  narrow  pulse. 

In  Section  7 we  considered  the  coupling  between  the  photon  field  and 

the  electronic  excitation.  There  we  concluded  [see  Eq.  (50)]  that  for  a 

given  peak  stimulated  emission  rate,  the  maximum  number  of  photons  generated 

in  the  cavity  is  proportional  to  the  lifetime,  therefore  one  would  like  long 

lifetime  centers.  Indeed,  in  the  case  of  Q-switching  we  found  that  long 

lifetime  centers  were  desirable.  However,  in  the  case  of  a two-center 

laser  the  lifetime  to  be  substituted  for  t in  Eq . (50)  is  an  effective 

lifetime  x , : 

a 1 
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From  the  considerations  of  Section  4 it  is  obvious  that  ^ t^,  i.e,, 
the  lifetime  of  the  fast  centers. 

In  Section  10  we  obtained  Eq.  (77)  which  gives  the  maximum  stimulated 
emission  rate.  The  lifetime  dependence  of  Eq.  (77)  results  from  Eq.  (74) 
which  can  be  written  in  the  form: 


where  this  effective  lifetime  is  given  by: 


1 


1 


(105) 


(106) 


Again  the  performance  is  dominated  by  the  fast  centers.  Having  exhausted 
effects  due  to  lifetime,  we  must  conclude  that  adding  slow  centers  has  no 
stong  effect  on  the  transient  spi)ce  operation. 

In  Figure  16,  we  use  numerical  integration  to  compare  the  transient 
spi)ce  for  a two-center  laser  having  lifetimes  differing  by  a factor  of  ten 
with  spi)ce  operations  for  individual  lasers  with  single  centers  of  either 
lifetime.  The  pea)c  intensities  for  the  three  are  comparable  and  no  apparent 
advantage  of  the  two-center  laser  is  shown.  One  might  suggest  that  some 
different  set  of  parameters  for  the  integration  might  demonstrate  enhanced 
intensity  for  the  two-center  laser.  However,  the  analytical  discussion 
offers  no  encouragement.  From  this  analysis  it  is  apparent  that  any  slow 
centers  that  are  present  can  store  energy  better  than  the  fast  centers  but 
in  transient  spilce  operation  merely  sit  on  the  sidelines.  If  two  types  of 
centers  are  present  in  a given  laser  material,  they  should  be  converted  to 
a single  lifetime  if  possible. 
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13.  CONCLUSIONS 


We  have  succeeded  in  obtaining  approximate  analytic  expressions 
describing  a wide  variety  of  laser  pulse  phenomena  that  can  be  handled  using 
the  laser  rate  equations.  Using  these  expressions  we  can  examine  these 
phenomena  without  recourse  to  the  catalogs  of  numerical  computer  data  that 
have  been  the  traditional  approach  to  laser  pulse  phenomena. 

In  Table  4 we  compare  various  pulse  modes  among  themselves  and  with 
cw  operation.  For  laser  materials  with  long  spontaneous  lifetimes,  the 
Q-switch  mode  is  to  be  preferred  since  it  gives  the  highest  peak  intensity. 
The  enhancement  of  the  peak  intensity  over  cw  is  approximately  equal  to  the 
ratio  of  the  spontaneous  lifetime  to  the  cavity  lifetime. 

As  the  spontaneous  lifetime  approaches  the  cavity  lifetime,  the  pre- 
ferred mode  is  cavity  dumping,  if  the  residual  cavity  losses  are  small 
enough.  In  this  case  the  enhancement  over  cw  is  proportional  to  the  ratio 
of  the  cavity  lifetime  with  zero  output  coupling  to  the  cavity  round-trip 
time . 

Transient  spike  operation  is  useful  in  roughly  the  same  range  as 
Q-switching.  However,  the  enhancement  factor  over  cw  operation  is  approx- 
imately equal  to  In  q^,  i.e.,  the  logarithm  of  the  number  of  modes  in  the 
luminescent  line  width.  Using  a typical  value  for  q^  gives  F % 22.  This 
is  low  relative  to  that  achievable  with  Q-switching.  Transient  spike 
operation  is  useful  if  one  is  interested  in  moderate  pulse  intensities  and 
in  dispensing  with  extra  intracavity  elements  or  if  one  is  principly  in- 
terested in  obtaining  narrow  pulses.  If  there  are  lasing  centers  of  dif- 
ferent lifetimes  present,  transient  spike  operation  is  dominated  by  the 
fast  centers  and  no  enhancement  of  peak  intensity  results  from  adding  long 
lifetime  centers  for  energy  storage. 
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SECTION  V 

SPECTROSCOPY  OF  A HIGH  CONCENTRATION  Al-PHOSPHATE  GLASS 


ABSTRACT 


An  Al-phosphate  glass  containing  Nd^'*’  In  concentrations  ranging  from 
20  21 

3 X 10  to  2.7  X 10  lons/cm  has  been  prepared  and  Investigated  spec- 
if 4 

troscoplcally . Different  ways  of  measuring  the  Fo/o  "*■  emission 

—20  ^ y/ £pq  2 

cross  sections  yield  values  between  1.96x10  and  2.7  x 10  cm  . At 
the  highest  concentration  the  decay  time  Is  50  ys,  while  the  radiative  life 
time  Is  estimated  to  be  > 446  ys.  Lasing  experiments  were  performed  by 
pumping  thin  platelets  of  glass  coaxially  with  a dye  laser.  Cross  sections 
losses  and  differential  efficiencies  are  derived  from  these  experiments. 

The  status  of  high  concentration  glasses  vs.  stoichiometric  Nd  compounds  Is 
viewed . 


PiOS 


SUiOMIOT 


1 . INTRODUCTION 


In  the  last  few  years,  there  has  been  considerable  activity  in  attempts 

to  make  small  lasers  which  Involve  stoichiometric  Nd  compounds  which  contain 

up  to  30  times  the  Nd  content  of  conventional  Nd  laser  materials  such  as 

Nd:YAG.  This  work  has  been  recently  reviewed  by  Chinn,  Hong,  and  Pierce^ 

2 

and  by  Danielmeyer.  The  most  promising  single  crystalline  materials  are 
LlNdP|^0^2  NdP^O^;^,  otherwise  known  as  meta  and  ultraphosphates  (the 

ultraphosphates  are  also  referred  to  in  the  literature  as  pentaphosphates ) . 

The  above  compounds  make  use  of  particular  crystalline  lattice  con- 
straints which  isolate  Nd  ions  so  that  typical  Nd-Nd  distances  are  in  the 
0.5  - 0.7  nm  range.  The  phosphate  materials  accomplish  this  by  formation 
of  polymer  chains  in  which  Nd  ions  are  separated  by  -O-P-0-  linkages. 

In  view  of  the  success  with  crystalline  stoichiometric  materials  the 

question  has  been  recently  raised  whether  development  of  high  concentration 

3 H 

glasses  is  another  possible  alternative.  ’ This  work  is  stimulated  as 

much  by  the  desirability  to  develop  miniaturized  Nd  lasers  as  by  the  use  of 

5 6 

glass  in  fiber  technology.  ’ One  possible  approach  is  to  form  glasses  of 
composition  close  to  the  meta  and  ultraphosphates.  The  Russian  workers 
have  been  quite  successful  in  obtaining  glasses  of  composition  close  to  the 
metaphosphate,  whereas  ultraphosphate  glasses  have  been  difficult  to 
control  and  are  hygroscopic.* 

Prior  to  the  advent  of  stoichiometric  materials,  several  groups,  pre- 
dominantly in  the  USSR  and  France,  have  been  active  in  preparing  various 

7 8 

low  concentration  Nd-doped  phosphate  glasses  for  large  laser  applications.  ’ 
The  limitations  which  emerged  from  this  work  were  due  to  concentration 
quenching  and  nonradlatlve  losses  attributable  to  OH  contamination.  The 
latter  effect  is  quite  analogous  to  our  results  on  liquid  hosts  for  Nd 
(namely  POCl^)  which  are  extremely  sensitive  to  water  content.^ 

A number  of  procedures  have  been  reported  to  reduce  the  water  content 
of  phosphate  glasses  by  as  much  as  two  orders  of  magnitude. These  in- 
clude bubbling  N^  through  the  liquid  during  founding,  addition  of  Cl  or  F 
as  getters,  a prolonged  founding  procedure,  and  evacuation  of  the  sample 
during  founding.  In  general,  however,  this  work  has  been  limited  to  glasses 
with  low  Nd  concentrations. 

In  totally  unrelated  work  at  GTE  Laboratories,  some  aluminum  phos- 
phate glasses  capable  of  accepting  rather  large  concentration  of  rare  earths 

• Private  communication,  J.  Gualtlerl,  Ft.  Monmouth,  New  Jersey. 


have  been  developed.  Normally  trlvalent  cations  such  as  aluminum  are  not 
considered  the  best  choice  for  laser  glasses  since  they  Invariably  lead  to  a 

O 

rather  broad  fluorescence  line  width.  In  this  investigation  we  were  not 
guided  by  the  choice  of  best  possible  glass  composition  but  rather  by  the 
availability  of  the  material  and  need  to  establish  a standard  of  reference 
for  further  work.  The  usefulness  of  such  a reference  Is  obvious  when  one 
considers  that  In  the  literature  of  phosphate  glasses  the  details  of  the 
spectroscopy  needed  to  fully  describe  a system  are  often  omitted. 

In  this  paper,  we  describe  some  spectroscopic  results  and  laser 
operation  with  high  Nd  concentration  aluminum  phosphate  glasses.  The 
results  are  encouraging  considering  that  no  special  attempts  were  made  to 
control  or  reduce  water  contamination . 


I 

I 

I 
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^ 2.  MATERIAL  PREPARATION 

f Neodymium-doped  aluminum  phosphate  glasses  have  been  prepared  by  heating 

mixtures  of  Al(OH),  and  Nd202  In  ~ 1 M H^PO^^  to  Increasing  temperatures  In 
order  to  remove  the  water  and  to  obtain  a homogenous  starting  batch.  The 
solid  mix  Is  fired  In  AI2O2  crucibles  at  1500°C  for  1 hour.  The  Initial 
composition  Is  chosen  to  correspond  to  the  metaphosphate  AlP^O^;  during 
firing,  P2O5  volatilization  shifts  the  composition  to  the  desirable  glass- 
forming region  In  the  vicinity  of  the  AlP^Og-AlPO;^  eutectic  which  Is  shown 
In  Figure  1.  After  pouring  Into  graphite  molds,  the  samples  are  annealed 
at  ~ 700°C.  Three  samples  of  aluminum  phosphate  glass  were  prepared  in  this 
manner  containing  31  mole  % (24  weight  %)  Al^O^,  and  their  Nd  content  Is 
listed  In  Table  1.^^ 


Figure  1.  The  Binary  System  Al^O,  • 3P,0c-  - AloO,  • PoOc 
(from  Ref.  12)  ^ 3 25  23  25 
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TABLE  1 


I 
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3.  EXPERIMENTAL  RESULTS 

3.1  EMISSION-ABSORPTION 

For  spectroscopic  purposes,  plates  of  5 x 5 mm  dimension  and  thickness 
ranging  from  0.025  cm  to  0.6  cm  were  cut  from  the  glass  buttons  and  polished 
on  both  sides.  Emission  spectra  were  taken  by  using  either  Xe  arc  or  dye 
laser  excitation  as  shown  on  Figure  2.  Decay  times  could  be  obtained  using 
the  same  apparatus.  In  this  case,  the  fluorescence  In  the  1.06  micron 
region  was  Isolated  by  an  Interference  filter  and  detected  by  a 
photomultiplier  (PMl). 

The  apparatus  could  further  be  modified  to  obtain  absorption  spectra 
at  low  temperature  which  are  needed  for  a correct  assignment  of  levels. 

A tungsten  filament  lamp  was  Imaged  on  the  sample  holder  which  In  turn  was 
Imaged  on  the  entrance  slit.  From  runs  made  with  and  without  a sample,  a 
transmission  curve  (not  corrected  for  reflection)  could  be  obtained. 

Flp-iire  3 gives  emission  spectra  at  room  temperature  of  the  three 
glass  sai..^  ..es  excited  by  monochromatic  radiation  from  the  dye  laser.  We  see 
that  there  Is  considerable  difference  between  samples  2 and  3 In  the 
resonance  region  (850  - 930  nm)  due  mostly  to  self-absorption.  The  shape 
of  the  laser  transition  line  (1030  - 1080  nm)  seems,  however,  hardly  af- 
fected by  the  Increase  In  concentration.  A spectrally  connected  emission 
intensity  (using  a standard  lamp)  for  sample  2 Is  shown  on  Figure  4. 

Figure  5 gives  the  absorbance  (=  a d log^^e;  a-absorptlon  coefficient, 

d thickness)  for  sample  2 taken  on  a Cary  14  spectrometer.  The  labeling 

of  the  states  Is  taken  from  Ref.  14.  In  making  spectroscopic  assignments 

of  the  various  lines  occurring  in  the  emission  spectrum,  a room  temperature 

absorption  spectrum  of  the  kind  shown  on  Figure  5 Is  of  limited  usefulness 

4 4 

because  the  transitions  ^5/0  ^3/2  r’^^olved.  In  order  to  obtain 

better  information  on  the  structure  of  the  ground  state,  low  temperature 
absorption  spectra  had  to  be  obtained. 

4 2 

Figure  6 gives  the  Ig/2  ^1/2  absorption  spectru'.i  taken  at  liquid 

nitrogen  temperature.  Since  the  upper  state  Is  unsolit  by  the  crystal 

J4 

field,  this  absorption  should  mirror  the  structure  of  the  Ig/2  ground 
state.  At  77°K,  It  does  not  seem  possible  to  observe  all  the  five  stark 
components  which  may  be  expected  to  be  present.  From  this  spectrum,  we 
can  only  Identify  the  components  Z^,  of  the  I9/2  manifold.  The 

levels  and  R^^  discussed  below  belong  to  the  ^^11/2  ^^3/2 

respectively . 


I 

I 


135 


Figure  2.  Schematic  Layout  of  Apparatus  for  Measuring  Lifetimes 
and  Time  Resolved  Spectroscopy 


Figure  4.  Spectrally  Corrected  Emission  Spectrum  of  Sample  2.  For  meaning  of  Shaded  Portions  see 
Section  3-3 
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The  splitting  of  the  ^2/2  tie  partially  obtained  from  a low 

temperature  spectrum  taken  In  the  region  around  875  nm.  This  Is  shown  on 
Figure  7-  Since  the  lines  are  still  quite  broad  at  this  temperature  due 
to  contribution  from  the  state,  the  ^3/2  splitting  shown  on  the  figure 
as  130  cm  ^ cannot  be  regarded  as  highly  accurate.  A better  value  will  be 
obtained  by  combining  these  data  with  low  temperature  emission  spectrum 
shown  on  Figure  8. 


It  Is  well  known  that  water  (or  OH)  contamination  of  the  host  material 
Introduces  nonradlatlve  quenching  of  the  luminescence.  In  order  to  obtain 
some  Information  on  the  water  content,  we  have  run  Infrhred  spectra  of  the 
samples  In  the  region  of  2.5  to  A microns.  An  example  of  such  a spectrum 
is  given  In  Figure  9.  The  band  at  3AOO  cm~^  marked  with  an  arrow  Is  at- 
tributed to  OH.  From  the  peak  transmission  of  82%  (not  corrected  for 
reflection),  we  obtain  an  absorption  constant  a (3^00  cm~^)  < 7-9  cm~^. 

This  value  will  be  compared  In  the  next  section  with  measurement  quoted 
in  the  literature. 


3.2  SPECTROSCOPIC  ASSIGNME.MTS 

From  the  data  of  Figures  3,  6 and  8,  we  can  arrive  at  a plausible 

assignment  of  some  six  lines  In  the  resonance  region.  Assignments  In  the 

laser  transition  region  are  more  difficult  because  the  lines  are  poorly 

resolved  and  we  have  no  Independent  information  on  the  structure  of  the 

ground  state.  There  Is,  however,  little  doubt  that  the  peak  occurring  at 

1052  nm  occurs  from  the  lower  ^2/2  Comparing  Figures  3 and  8 

we  see  that  this  peak  remains  largely  unaffected  by  temperatures,  whereas 

the  short  wavelength  shoulder  of  the  emission  band  sharpens  considerably, 

leaving  a small  residual  peak  (transition  from  R^)  at  77°K  and  located 

at  1036  nm.  The  energy  difference  between  these  peaks  corresponds  to 

147  cm~^  and  hence  Is  within  10  cm~^  of  the  R^  - splitting  determined 

from  the  resonance  transition  data  (Figure  6).  The  summary  of  the  asslgn- 

A 

ments  Is  given  on  Figure  10  and  In  Table  2.  The  extent  of  the  ^11/2 
state  on  the  long  wavelengths  side  (Yg  level)  can  only  be  given  approx- 
imately because  of  the  tailing  off  response  of  the  SI  photomultiplier. 

3.3  EMISSION  CROSS  SECTION 

From  the  soectroscoplc  data,  we  can  determine  the  emission  cross 

4 4 

section  for  the  transition  ^2/2  ^9/2'  There  are  several  ways  of 

handling  this  problem,  resulting  In  somewhat  different  expressions  for 
the  cross  section.  In  each  case  the  basic  Idea  Is  to  compare  the 
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Figure 
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Figure  9.  Transmittance  Spectrum  of  Sample  1 In  the  Infrared 
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fluorescence  Intensities  of  the  laser  ( F,/o  -►  In/o)  transition  and  the 
4 if 9/2 

resonance  ( F2/2  ^11/2^  transition  for  which  the  absorption  cross 

section  can  be  measured  directly. 


TABLE  2 


WAVELENGTHS  AND  ASSIGNMENTS  OF  TRANSITIONS 
X(nm) 


865 

877 

873 

885 

892 

905 


^•=3/2 


9/2 


1036 

1052 


^3/2  - *1 


1/2 


A straightforward  application  of  the  Ladenburg-Fuchtbauer  equation  for  the 
integrated  absorption  of  a spectral  line  leads  to  the  following  expression 
for  the  peak  emission  cross  sectlon^^ 


The  subscripts  n and  m refer  to  the  initial  and  final  states  in  the 

emission  process.  The  subscripts  1 and  J refer  to  the  initial  and  final 

states  of  the  absorption  process  characterized  by  the  peak  absorption  cross 
A 

section  a . The  wavelengths  A designate  the  peak  of  the  emission  and 
absorption  lines;  N's  are  the  populations  of  the  states  and  g's  are  their 
degeneracies.  The  quantities  I refer  to  the  peak  intensities  of  the 
emission  (watts  per  unit  wavelength  interval)  for  the  n-^m  and  J-^1  lines. 

Eq.  (1)  Involves  no  approximations  other  than  the  equality  of  local 

14 

field  corrections  for  the  various  states  and  is  well  suited  for  the 
calculation  of  cross  sections  if  the  n-*-m  and  1->J  transition  can  be  unam- 
biguously assigned.  Such  is  the  case  of  Nd  in  YAG  to  which  Eq.  (1)  was 
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applied  In  Ref.  13.  Use  of  Eq.  (1)  for  a substance  like  glass  Is  not  en- 
tirely Justified  because  the  lines  connecting  Individual  Stark  components 
of  the  F-,/-,,  In/-,  and  /n  states  are  not  well  resolved.  It  Is,  In 


3/2’  "9/2 


11/2 


,15 


principle,  possible  to  follow  the  procedure  of  Mann  and  DeShazer  and 
decompose  the  broadened  spectrum,  but  such  procedure  Is  both  tedious  and 
not  necessarily  unique.  With  this  proviso  It  Is  nevertheless  Instructive 
to  apply  It  to  the  glass  case.  From  the  spectroscopic  assignment  It  ap- 
pears reasonable  that  the  peaks  at  877  nm  and  1052  nm  correspond  to  the 

If  we  assume  that  other  transitions  con- 


R, 


and  R^  -♦  transitions. 


tribute  negligibly  to  the  peak  Intensities  at  these  two  wavelengths  (a  fact 
schematically  Illustrated  by  the  shaded  portions  on  Figure  4)  then  we  can 
use  the  ratio  of  these  two  Intensities  directly  In  Eq . (1).  We  further 
need  the  absorption  cross  section  for  sample  2 corresponding  to  the  R^ 

transition  (877  nm).  The  peak  absorption  constant  Is  measured  to  be 
2.3  cm”^  for  this  line.  The  peak  cross  section  Is  then  given  by: 


peak 

nTzTT 


(2) 


where  the  population  of  the  state  Z^  Is: 


N(Z^) 


gqNo® 


-E^/kT 


'S^e 


-E^/kT 


(3) 


,3+ 


Nq  Is  the  concentration  of  Nd  per  cc  and  g^^  Is  the  degeneracy  of  the 
state.  Using  the  three  Z states  from  Figure  10  and  assigning  to  each  a 
Kramers  degeneracy  of  2,  we  obtain  for  the  partition  function  a value  of 
3.59  and  for  the  peak  absorption  cross  section: 


a^CZi-vRi)  = 0.5  X 10  cm^ 

.Since  both  emission  and  absorption  transitions  involve  the  same  R, 


level : 


From  these  values  we  obtain  from  Eq.  (1); 


(5) 


a‘"(Rl--yi) 


1.96  X lO"^'^  cm^ 


(6) 
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There  exists  another  way  of  calculating  the  emission  cross  section, 

widely  used  in  the  glass  llterature^^  and  first  formulated  by  Sarkies, 

17 

Sandoe  and  Parke.  In  this  formulation  the  individual  and  unresolved 
Stark  components  of  a state  are  all  lumped  together  and  the  cross  sections 
(both  for  emission  and  absorption)  are  averaged  for  the  whole  state.  The 
result  is; 


5^ 


(7) 


where  the  subscripts  1 and  r refer  to  the  laser  and  resonance  regions;  the 
factor  5/2  is  the  ratio  of  degeneracies  of  the  ground  and  metastable  level 
( ^^eff  defined  by; 

/^.(DdX  (8) 

A.  \ _ 

eff  ~ I^(peak)  ‘ 

The  Integrated  absorption  cross  section  appearing  In  Eq . (7)  Is  obtained 
from  an  absorption  spectrum  similar  to  that  of  Figure  5 but  taken  on  a 
sample  of  0.188  cm.  We  obtain: 


0.6H  X 10"^^  cm^  . (9) 


With  = 28.8  nm  the  emission  cross  section  is: 

The  discrepancy  between  the  two  values  of  the  emission  cross  section 
Indicates  the  sensitivity  of  the  result  to  the  method  of  calculation,  the 
assumptions  that  go  Into  the  use  of  Eqs.  (1)  and  (7),  and  Indeed  the 
definition  of  the  quantity  which  Is  often  loosely  called  the  "emission 
cross  section." 

1 B 

It  has  been  pointed  out  recently  that  the  averaging  which  goes  into 
the  derivation  of  Eq.  (7)  may  be  seriously  at  fault  If  the  transitions 
between  Individual  Stark  components  of  Initial  and  final  states  are  not 
equal  and  If  their  populations  differ  significantly.  Strictly  speaking, 

Eq.  (7)  Is  only  valid  when  the  states  are  degenerate.  It  Is  thus  difficult 
to  pinpoint  exactly  the  reason  for  the  discrepancy,  but  the  results  Illus- 
trate the  pitfalls  of  applying  Indiscriminately  expressions  taken  from  the 
literature.  A possible  way  of  obtaining  more  reliable  values  of  emission 
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cross  section  has  been  indicated  In  Ref.  l8  but  this  requires  data  taken  at 
liquid  helium  temperature. 

A comment  should  be  made  about  the  role  of  self-absorption  In  the 
determination  of  the  cross  section.  Examining  the  spectra  of  Figure  3 It 
Is  obvious  that  the  spectrum  of  sample  3 Is  heavily  distorted  by  self- 
absorption In  the  resonance  region.  Samples  1 and  2,  however,  show  almost 
Identical  spectra  — although  they  differ  by  a factor  of  2.8  In  Nd  concen- 
tration. We  have  also  calculated  the  self-absorption  correction  for  sample 
2 and  concluded  that  It  may  account  for  an  error  In  the  resonance  spectrum 
of  less  than  ten  percent.  Neglecting  self-absorption  can,  therefore,  only 
result  in  a slight  overestimate  of  the  cross  sections. 


3.^  LIFETIMES 

The  luminescence  decay  times  were  measured  for  the  three  samples  by 
using  short  pulse  (dye  laser)  excitation  at  58^  nm  and  the  boxcar  averager. 
For  the  most  concentrated  sample  3,  a slight  deviation  from  exponential 
decay  Is  observed  In  the  early  stages  of  decay  as  shown  on  Figure  11.  Such 
a deviation  is  common  In  Nd  doped  glasses  and  is  a consequence  of  the 
statistical  distribution  of  Ion  sites. 

Heat  treatment  of  glass  3 at  750°C  for  2 hrs  In  vacuum  or  In  oxygen  did 
not  produce  any  measurable  change  In  the  lifetime.  To  the  degree  that  the 
decay  time  of  the  most  concentrated  sample  is  due  to  OH  contamination,  much 
longer  annealing  times  and  possibly  higher  temperatures  may  be  required  to 
reduce  this  contamination. 

The  dependence  of  decay  time  on  concentration  Is  presented  on  Figure 
12  (curve  marked  Al-Phosphate ) . Unfortunately,  we  are  not  certain  of  the 
limiting  value  of  the  decay  time  at  low  Nd  concentration  since  low  doped 
samples  were  not  available.  In  the  same  figure  we  have  also  entered  data 
on  other  materials  which  we  shall  discuss  In  the  next  section. 

i| 

The  radiative  lifetime  of  the  state  ^3/2  obtained  In  terms  of 

the  cross  section  for  a luminescent  transition  connecting  this  state  to  a 
lower  state  (say  the  branching  ratio  y for  this  transition 
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Figure  11.  Decay  of  Fluorescence  of  Sample  3 at  Room  Temperature 
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where  we  have  taken  n = 1.55.  ' An  accurate  determination  of  the  branching 

ratios  require  measurements  of  emission  intensity  for  transitions  terminating 

on  all  members  of  the  quartet  I manifold.  This  has  not  been  carried  out  for 
4 4 

the  states  ^13/2  ^15/2  ‘ Assuming,  however,  reasonable  values  taken 

from  the  literature , we  obtain: 

^S/2  ^ 0-385 

^11/2  = 0.514 

^13/2  i 0.1 

^15/2  “ ° • (12) 

From  Eq . (11)  and  these  values  we  obtain: 
j 446  ys 

I 6l4  ys  * (13) 

depending  whether  the  larger  or  smaller  of  the  two  cross  sections  were  taken. 
These  two  values  of  the  calculated  radiative  lifetime  could  be  checked 
against  measurements  of  decay  time  taken  on  very  dilute  samples.  On  the 
basis  of  available  data  (see  Figure  I6),  it  is  not  possible  to  say  which  is 
closer  to  the  true  value. 
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For  the  purpose  of  comparison,  we  present  In  Table  3 some  basic  data  on 
a few  representative  phosphate  glasses.  We  see  from  this  comparison  that 
parameters  of  our  aluminum  phosphate  glass  fall  well  within  the  range  of 
those  measured  on  other  preparations. 

To  our  knowledge,  only  the  two  last  entries  (Ll-MP  and  the  Al-phosphate ) 

have  been  Investigated  at  Nd  concentrations  substantially  exceeding 
21  3 

10  lons/cm-^,  and  therefore  comparable  to  the  stoichiometric  Nd  compounds. 
Prom  Table  3 we  see  that  these  two  glasses  are  roughly  equivalent. 

TABLE  3 

COMPARISON  OP  SOME  PHOSPHATE  GLASSES 


LHG-5* 

P107* 

> EV1* 

L-41* 

Li-MP*» 

Al- 

Phosphate 

^^eff 

25.5 

26.5 

21.7 

32.5 

31.0 

28.8 

Xj  (nm) 

1054 

1054 

1054 

1056 

1055 

1052 

’’rad 

351 

322 

314 

564 

330 

446 

(10-20  ^^2, 

3.9 

3.9 

4.7 

1.9 

3.2 

2.7 

* Ref  16 

**  Li-Metaphosphate  Pef  4 


Another  perspective  on  high  concentration  phosphate  glasses  Is  obtained 

from  Figure  12.  From  It  we  see  that  both  the  aluminum  phosphate  and  Ll-MP 

19 

fall  between  the  concentration  quenching  curves  of  Nd  ultraphosphate  (UP) 
and  the  French  phosphate  glasses  reported  In  Ref.  8.  The  data  on  tellurite 
glasses  were  Included  In  the  figure  since  they  were  recently  reported^  and 
suggested  as  good  candidates  for  mlnllasers. 

The  quenching  characteristics  of  the  Al-phosphate  glass  are  essentially 
similar  to  the  Ll-MP  glass.  Since  very  dilute  samples  were  not  available 
at  this  time,  we  do  not  know  to  what  extent  the  decay  time  would  approach  the 
calculated  1)46  ps  radiative  lifetime  of  the  ^^2/2  Should  there  be  a 

substantial  discrepancy,  we  would  have  to  conclude  that  OH  quenching  occurs 
In  our  samples.  The  value  of  the  Infrared  absorption  constant  (8  cm”"^), 
derived  from  Figure  9,  compares  favorably  with  values  achieved  by  Soviet 
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Density  of  Nd  (cm 


Figure  12.  Concentration  Quenching  of  Some  Nd  Containing  Materials: 

• - present  A1  phosphate  glass;  X - Li  metaphosphate  glass 
from  Refs.  3 and  4;  "UP"  - Nd  ultraphosphate  from  Ref.  19; 
"French"  - phosphate  glasses  from  Ref.  8;  "Tellurite"  - glass 
from  Ref.  6 


154 


workers  who  used  halogenlzatlon , bubbling  and  long  founding  times,  since  none 

of  these  procedures  were  attempted.  It  Is,  however,  well  above  an  absorption 

of  1 cm~^  at  3-^5  microns  achievable  by  extra  care.^'^  This  leads  us  to 

believe  that  the  Al-phosphate  glasses  do  stand  a chance  of  considerable 

Improvement  In  that  respect.  We  find  It  Is  somewhat  surprising  that  Ll-MP 

has  much  larger  bandwidth  than  Al-phosphate,  since  according  to  the  work  of 

0 

Deutschbeln,  et  al..  It  appears  that  monovalent  cations  produce  the  narrowest 
lines,  whereas  trlvalent  cations  (specifically  Al)  produce  the  largest. 

For  practical  applications  the  chemical  stability  and  resistance  to 
moisture  can  be  of  major  significance.  In  that  respect  the  aluminum  phosphate 
glass  appears  to  be  excellent.  Over  a period  of  months  we  have  not  observed 
any  deterioration  of  the  polished  surfaces  under  ordinary  laboratory 
conditions . 

In  summary  our  results  Indicate  that  high  concentration  glasses  should 

not  be  neglected  In  the  search  for  iiew  Nd  hosts  since  concentration  quenching 

In  these  materials  plays  surprisingly  less  of  a role  than  generally 
19 

anticipated.  Finally,  when  the  technology  of  preparation  of  this  glass 
Improves,  the  possibility  of  making  large  melts  will  offer  a distinct 
advantage  over  one-at-a-tlme  crystal  growth. 

The  authors  are  grateful  to  Dr.  L.  Rlseberg  for  making  the  Nd  glass 
samples  available.  Many  helpful  clarifications  on  the  chemistry  aspects 
are  due  to  Dr.  B.  McCollum.  The  absorption  spectra  were  kindly  supplied 
by  Dr.  P.  Cukor.  Dr.  E.J.  Johnson  supplied  a computer  program  for  the 
correction  of  the  spectra. 
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